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ABSmVJT 


A theoretical study by J, E. Chance and J. M. Cantu is included in 
the Appendix of this paper on the vegetation model of G. H. Suits 
The experimental bidirectional reflectance of cotton is presented and 
compared to the Suits model. Some wheat reflectance data is presented 
for a Mexican dwarf wheat. The general results are that the exchange of 
source position and detector position gives the same reflectance measure- 
ment if the irradiance is purely specular. This agrees with Suits. The 
reflectance versus sun angle and reflectance versus detector angle do 
not agree with the Suits predictions. There is qualitative agreement 
between the Suits model and reflectance versus wavelength, but quantitative 
agreement has not been observed. Reflectance of a vegetation canopy with 
detector azimuth shows a change of 10 to 40% for even sun angles near 
zenith, so it seems advisable to include azimuthal angles into models of 
vegetation. 


( 1 ) 


G. H. Suits, Remote Sensing of Envirorment 2,117 (1972), 
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FINAL REPORT 


INTRODUCTION 

The task of Identifying varieties of agricultural crops from satellite 
altitude has proven to be a very challenging problem. An empirical appro- 
ach using statistical techniques has been used predominantly, but in 
recent years a physical model approach which includes both models of the 
atmosphere and of the plants themselves has been used to try to estab- 
lish cause-effect relationships. 

The number of temporal variables involved in the problem renders the 

empirical approach useless if a training field is used to categorize crops 

( 1 ) 

over a period as long as 1 or 2 hours as was sho^ra by Malila, et . al . 

At the same time these time-dependent parameters make physical models so 
complicated that they rapidly become formidable. 

The Suits model was an attempt to use only the predominant physical 
parameters of a plant canopy that affect the reflectance. This report is 
a study of the properties of the Suits model on cotton and, to a limited 
gxtent, wheat to detemaine if tUe dominant characteristics have, indeed, 
been modeled. The principal investigator aroused the interest of a pro- 
fessor from the Mathematics Department at Pan American University, Dr. 

J. E. Chance, and a Master's Degree Candidate, Mr. Juan Manuel Cantu, to 
further study the mathematical properties of the Suits model. As far as 
was practicable the same questions that were asked of the model were asked 
of the field measurements with comparisons being made in this report. 

The main factors that are of interest are the folloTrlng: 

(a) Will the reflectance measurement be affected by interchange 
of source and detector? 

(b) In what regime of source position and observer position is 
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the bidirectional reflectance function the sjost sensitive? The 
leai't sensitive? 

(c) Does the top surface reflection of leaves, which in many plants 
is admittedly non-Lamhertian, as observed in single-leaf experi- 
ments of Breece and Holmes account for the non-Lambertian 
properties of vegetative canopies. 

(d) In general characteristics, how is the qualitative agreement 
between the model and field measurements? Is the agreement 
about the same for low leaf area index (L.A.I.) crops and for 
high L.A.I. crops? 

(e) Is the azimuthal variation negligible as implied by Suits' 
assumption of azimuthal S3unmetry? 

Most of the data were collected on cotton, sampled from faro different 
fields where the L.A.I. was 5.3 and 5.5. The wheat was a Mexican dwarf 
variety gro;m in Eagle Pass, Texas. The variety was Penjamo and is nearly 
identical to the high-yield varieties introduced 10 years ago into India 
and Pakistan and now called Kalyan, Kalyan-sona, and Sona-lika. The L.A.I 
of the wheat reported here was 2.1 and was in a mature green stage mth 

well-developed heads. 


THEORY OF MEASUREMENTS 

iVl 

The reflectance is defined as p ~ , where or ]e^| are the 

magnitudes of the electric field vectors of the reflected and incident 
radiation, respectively. This can easily be T<rritten as the ratxo of the 
total time-averaged "outward" moving radiant power relative to the surface 
being considered. In a laboratory spectrometer single leaf reflectance is 


normally measured by irradiating the leaf with white light from a diffuse 
source and measuring the light reflecting in a particular direction. In 
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order to apply the definition of the reflectance, it appears to this inves- 
tigator that one should specify the type of source and type of detector as! 
"hemispherical - hemispherical reflectance" or "hemispherical - angular 
reflectance", etc. This would also dictate that the reflectance function 
should be defined as 

radiance 
P = irradiance 

as done by Breece and Holmes . 

In the case of field studies, the source of irradiance is a combination 
specular from direct sunlight and diffuse from skylight. In any field 
studies one should separately measure both components. 

The reflectance of the canopy as determined in this report is the 
bi-angular or angular-angular reflectance. The ratio of specular to diffuse 
solar irradiance varied from -95% to a worst case of 67%. If one assumes 
pure specular flux from the sun, then the reflectance can be determined as 

follows: / 

(ca nopy) f P ^ (sun) 

Pc ” ~A%et^ cos0Det / ^ ^°®^un 

where the crop is viewed xrf.th a detector looking through a solid angle 

at a projected area AA. cosSj^gj- 

and 

A^sunAA ^°^®sun 

is measured irradiance of the 

sun falling on an area AA. 

P^ is the power in watts 

arriving at the detector in 

the wavelength interval centered on X._ This assumes a completely black sky 
and a sun subtending a solid angle AQ^^^ at the position of the observer. 
Since for field measurements, a technique must be utilized that xntegrates 
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the. total iT 5 xadi,ance . of direct sunlight plus diffuse light. This can be 
accomplished in the following manner; Use a diffusely reflectxng panel 
having a kno^ra reflectance, (panel), placed in the same position as the 
canopy whose reflectance is to be determined. One then obtains for the 

panel 

(panel) = AA oos0 

where AA is the area seen by the detector. This then yields an integrated 

value for the total solar irradiance; the canopy reflectance becomes 

(canopy) / (panel) ^ 

p (canopy) AAjj^^ cos^gy/ pi (panel) ASijjg^. AAj^t°°®®Uet 

Because of the experimental difficulty of having a reflectance standard large 
enough to be observed at the same position as the agricultural canopy, a 
horizontal panel is placed about 1 meter below the detector head and observed 
from a vertical position of the head. This requires that in the denominator 


Sx (sun) 


Det 


AJ^sun ^ ^°s®sun 


®Det “ equation thus yields 

(canopy) 


Pj^ (panel) 

Mjjgt pFlpanel) 


p (canopy) = 


iCcanopy) / 


pi (panel) 


COS0 


'Uet 


A / ^et 

One thus has a measuring technique that does not involve the geometry of the 
detector. 


experimental METHODS: 

Reflectance measurements taken in the field are made from the top of 
portable construction platform using a wedge-filter type radxometer 
(Isco Model SR), equipped mth a 1.8 meter fiber optics probe. The field of 
view to half maximum is 13° and full field is 19°. The spectral band 
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width is 15 nanometers (nm) in the visible and 30 nin. in the infrared with 
a range of 380 nm to 1550 nm. 

Spectral intensity was measured in the field by viewing the canopy 
from 4,9 m above the top of the plants at some azimuth angle measured from 
0° as magnetic north and a polar angle measured from 0° straight down 
along a plum1:i line. The radiometer readings were combined VTith the absolute 
reflectance of the reflectance panel to give the absolute bidirectional 
reflectance function for the canopy as detailed in the previous section. 

The reflectance panel was spray painted with at least 3 layers of 
Eastman Kodak Tilhite Reflectance Coating #6080 made of barium sulfate and 
having a reflectance >97%. The absolute reflectance was determined perio- 
dically from the United States Department of Agriculture, Agricultural 
Research Service (USDA, ARS) at Weslaco, Texas, standard vitrolite sample 
using a Beckman DK-ZA*'^^. 

The variations in radiance from the canopy as a function of azimuthal 
and of polar angle were recorded on an X— Y recorder as folloX"7st A drafting 
TTiarhine was disassembled and modified to allow the radiometer fiber optics 
probe to be clamped at any polar angle. A 0,1% precision potentiometer 
X- 73 S placed above the rotating head as shoxm in Fig. 1 to give a voltage 
read-out proportional to the angle displayed on the X-axis of the recorder 
while simultaneously the Y-axis is the radiometer output. This 180° 
azimuth is scanned in about 1 minute and reversed to give some indication 
of the noise and transient fluctuations, 

RESULTS And DISCUSSION: 

(a) The first result is the effect of interchanging the source and 
detector positions. The altitude or polar angles were measured from zenith 
and the azimuth angles are magnetic compass measurements. 



TOP VIEW 



SIDE VIEl^ 


rTTW-r* 

mnn-( 



Figure 1« Azimuthal scan, device made from a drafting machine. 

The position sensing mechanism is a 1 turn lOOOfi helipot 
connected to a 2.2 volt battery. 




6 

TABLE 1. Effect of exchanging sun and detector positions. 

Cotton measurements are shovm for run prefix C and 
wheat is prefix W. 


RUN 

DATE 

SUN POLAR 
& AZIMUTH 

DETECTOR POLAR 
& AZIMUTH 

WAVELENGTH 

(nm) 

REFLECTANCE 

C50 

7/12/75 

49°, 270° 

45°, 90° 

500 

0.020+0.002 

C7 

7/21/75 

55°, 90° 

45°, 270° 

500 

0.020+0.002 

C51 

7/12/75 

50°, 270° 

45°, 90° 

850 

0.22+0.06 

C8 

7/21/75 

54°, 90° 

45°, 270° 

850 

0.27i0.06 

W39 

4/19/75 

35°, 270° 

30°, 90° 

850 

0.37+.015 

W61 

4/19/75 

27°, 90° 

30°, 270° 

850 

0.30+.014 


The indicated errors are the combined effects of fluctuations on the radio 
meter signal from the cotton crop and from the reflectance panel and the 
estimated error in the absolute reflectance of the standard panel. The spec- 
ular flux was 78% of the total at 850 nm and 75% at 500 nm on 7/21/75 at a 
solar polar angle of 45°; on 7/12/75 the specular flux was 89% at 850 nm and 
83% at 500 nm at a solar polar angle of 11°. The fact that the solar polar 
angles were so much different for these measurements makes it impossible to 
compare sky conditions other than to say that both days were quite clear 
(based on similar data reported by Jones & Gondit^^^. Total solar irra- 
diances for the two days in question were 233 watts/cm^ for 7/12/75 and 225 
on 7/21/75 at 850 nm; these differ by only 4%. This is justification for 
considering the sky conditions essentially the same. The 500 nm data has 
small error bars and shoxjs exchange s 3 nnmetry; the 850 nm data has large 
error bars and shows differences smaller than error limits. We thus con- 
clude that the exchange of source and observer gives reflectance values that 

are essentially the same for cotton. 

The exchange of source and observer for xjheat is shoxTO. There is not 
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perfect exchange because the sun angle was 35° and 27° for the cases shown. 
The specular flux was variable because of the presence of high cirrus 
clouds; it varied from 61% to 68% specular. Because of these factors, the 
exchange property for wheat cannot be verified, 

(b) The variation of reflectance of cotton with sun angle is shown in 
Fig, 2 and Fig, 3 for two angles of view at 500 nm and 850 nm, respectively. 
The general trends are toward moderate variations until the sun angle is 
more than 50° from vertical; the variation in reflectance is more pronounced 
in the near infrared compared to 500 nm. The data indicate that a minimum 
exists in the 850 nm data for a sun angle of about 40°-50^ in cotton. The 
Suits model prediction is shown in Fig, 7 of the report of Chance and Cantu 
in the Appendix of this paper. The saddle-point behavior predicted by 
Suits shows a "plateau" and decreasing reflectance versus sun angle for a 
view angle of <35°, In the observer angle range >50°, there is an increasing 
reflectance versus sun angle; in the intermediate range the reflectance is 
nearly flat. 

The only qualitative comparison one can make betoeen the Suits pre- 
diction and experimental values is that there is a change in the general 
behavior of the reflectance in the vicinity of observer angle 40°, Since 
the cotton canopy is nearly uniform, it is clear that the effect is not a 
result of a layer structure at the tops of the plants reflecting most of the 
light at extreme solar angles and extreme observer angles. If one accepts 
the exchange symmetry, it follows that the isoreflectance surface must have 
diagonal symmetry and can have a maximum, a minimum, a saddle— point or 
completely flat. This all assumes the relative insensitivity to the azimuth 
angle between source and observer. We proceed to now look at reflectance 
versus observer angle at fixed sun angle, i.e., slicing the isoreflectance 
surface in the direction of fixed sun angle. 




Figure 2. The reflectance of cotton at 500 nm for various sun polar angles. 

The detector polar angle was at^0° and x 30° looking at an azimuth 
angle 180° from the sun azimuth. 


REFLECTANCE VS" SUN ANSLE 



SUN ANGLE 


Figure 3, The reflectance of cotton at 850 nm for various sun polar angles 
The detector polar angle was at O*’ and x 30° looking at an 
azimuth angle 180° from the sun azimuth. 
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(c) The observer angle influence on the reflectance is observed in 
Figs. 4 and 5. The reflectance is shown here to be very sensitive to the 
polar angle of view at 850 nm, while the variation at 500 nm is only slight. 
The vertical plane through which the detector was scanned contained the sun 
in every case (the effect of choosing other planes of observation can be 
found from the azimuthal scans). 

The combination of reflectance versus sun angle and reflectance versus 
observer angle makes possible an eKperimental plot of isoreflectance curves 
for purposes of comparing these results with the Suits model predictions. 

The general features of this plot for cotton do not agree with the saddle- 
point behavior predicted by the Suits model. The slices at constant sun 
angle and variable detector angle show maxima as seen in Figs. 4 and 5. 

This may be attributable to either top surface specular reflectance from 
the leaves or to heliotropic movement of the leaves. More experimental 
evidence will be presented to support the first effect; the second effect 
seems to explain the maximum that occurs in the data presented in Fig. 4 
and Fig, 5, The maximum reflectance occurs when the observer has his back 
to the sun looking at about -5° in the Ik and at -30° in the visible. 

This general behavior of the reflectance might be expected for a collection 

of leaves oriented toward the sun. 

Shown in Figs. 6 and 7 is the reflectance of wheat at constant sun 
angle and variable detector polar angle. The negative sun angle indicates 
a morning sun. Note that there is no sign of row effects because the rows 
were only 0.16 meters apart. Since the L.A.I. was only 2,1, it seems that 
the data may have been greatly affected by the light, sandy-colored soil 
background. Looking vertically downward, the soil was readily Trisible; 
whereas at an angle of 30° or more only vegetation could be seen. The bare 
soil reflectance at 500 nm was 0.13 and at 850 nm it was 0.31* 

The effect observed by Breece and Holmes on single-leaf non-Lambertian 


SCAN #27 


7/12/75 


a 

54 n 





Figure 4* 


a 


PSLAR ANSLE 

Cotton reflectance measurements at 500 nm at a sun polar angle of , i 

31^ from zenith in the west. The detector polar angle is posxtxve when lookxng 

west in the vertical plane containing the sun. 


*~S* Q IS 

Pi3LAR ANSLH 


Figure 5. Cottoii reflectance at 850 nni tTith a sun polar angle of 32° front 
zenith in the west. The positive direction is toward the west. 




REFLECTANCE 

»2QQ .sea 


H-7P t SUN AN6LE=-40 850NM 



Figure 6. T'Jheat reflectance at 850 nm for a sun angle of 40° fror* zenitli in 

the eastern sky. The positive polar angle indicated is the detector 
polar angle directed toward the north. The irrigation network in the 
field would not ailov/ east-west polar scans. 


SOP 5 SUN AN6LE==-38> SOONH 



Figure 7. Tflheat reflectance at 500 nin for a sun angle of 38° from zenith in 
the eastern sky. 
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reflectance was expected from the top surface of the canopy observed at 
extreme sun angles. Observer polar angle versus reflectance is sho^m in 
Figs. 8 and 9, The negative angles for the sun were morning positions; the 
negative angles for the observer polar angle were with the observer looking 
westward in the plane of the sun. The reflectance showed a high value when 
the observer was looldng sunward in both morning and afternoon at 500 lun. 

The effect was not seen at 850 nm. This agrees mth the observations of 
Breece and Holmes on soybean leaves with a source polar angle of 60® 

(shown in their Fig. 6) . Their observation was that there was a. very strong 
tendency for specular reflectance in the visible wavelength ranges and the 
IR looked nearly Lambertian. One interesting feature of the 850 nm data 
was the marked effect of row structure at extreme sun angles compared to 
only a small effect at 500 nm (See Figs. 10 and 11). Compared to Figs. 8 
and 9, the row effects are much smaller, while the sun has moved only 12®, 
These results can be accounted for by the following discussion. Cotton has 
very large, nearly flat leaves. These leaves are nearly symmetrically dis- 
tributed about the central stalk of the plant as we determined by stripping 
5 randomly-selected plants. The average leaf slope was 32°. The helio- 
tropic effect was observed with time-lapse movie photography, but the mag- 
nitude of the change in leaf position was only a few degrees in the leaves 
in the upper part of the plant. The asymmetrical reflectance versus observer 
angle curves were probably due to non-Lamb ertian leaf reflectance as observed 
by Breece and Holmes on corn and soybeans. The row structure seen in the 
IR and not in the visible was most lilcely due to the fact that the IR pene- 
tration into the canopy is about 7 leaves deep and in the visible it is only 
about 2 leaves deep^^^. Therefore, the area between the rows which had a 
low L.A.I. of 2 to 5 would contrast with that in the rows of 10 or more if 
the light was penetrating to 7 leaves deep as in the infrared. To the human 
eye the canopy was nearly uniform. 




Figure 8. Cotton reflectance at an extreme sun angle of -79° at 500 nm. 

The scatter was not noisej data were repeatable and showed 
sinusoidal trends presumably from shadows between rows of plants 
The positive angles sho\?n are looking eastward toward the sun. 


SCAN #6 


7 / 31/75 



PGLAR ANSLE 


Figure 9. Cotton field reflectance at 850 nm. The spikes are row effects 
and were smooth, nearly sinusoidal curves in the original analog 
data. The sun angle was at -79° and the detector was looking 
eastward at positive polar angles. 





ABSOLUTE REFLECTANCE 


SCAN 7/31/75 



-»0 -•« -J6 O IS 30 4S 0O 


POLAR ANGLE 

Figure 11. Cotton reflectance at 850 nm for a sun angle of -68°. 

Positive direction for the detector is eastward toward 
the sun. 
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(d) The qualitative agreement between the data gathered in the field 
and the Suits model is discussed beginning on page 39 of the report in the 
Appendix of this paper. 

(e) The azimuthal variations in reflectance are shoitvm for various sun 
angles in Figs. 12 and 13 for cotton and in Figs. 14 and 15 for wheat. The 
only features are row structures that show in the 850 nm data for cotton. 
Azimuthal positions of the detector yield different values of the measured 
reflectance by as much as 25 % (Fig . 14) , so one is forced to admit that 
azimuthal angles are important when modeling vegetation. 

CONCLUSIONS: 

1. Exchange symmetry predicted by the Suits model for interchange of 
source and detector is generally verified for pure specular irradiance, 

2. The bidirectional reflectance function for cotton shows large 
variation in the IR at extreme sun angles (i.e., >50° from zenith) and 
moderate variation in the visible. 

3. The general behavior of the bidirectional reflectance function 
with sun angle and observer angle predicted by the Suits model does not 
agree with the observed experimental data for cotton. 

4. Azimuthal variation is of the order of 10% — 40° for both cotton 
and wheat. 

5. The wheat reflectance data shows very few distinct trends. No 
effects of row structure are apparent. Insufficient data exist to make 
generalisations about trends of the reflectance with sun angle and observer 
angle. 



ABSOLUTE REFLECTANCE 


SCAN #37 


7 / 12/75 



■ Figure 12. Cotton reflectance at 830 mn for a detector polar angle of dS" and a 
variable azimuth ahotm as the direction faced by the_detector, i.e., 
from south through west to north. The sun was at 42 from z-enlth 
in the west. 
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Figure 13. Cotton reflectance at 500 nm for detector polar angle of 45°, variable 
azimuth. The sun xiras at 42° from zenith in the x^est. 
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cluded in Appendix 1. Also to be found in /{jpendix 1 
are Siaits’ parameters for cotton and wheat. 
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1. Introduction 

Many of the problems relating to the growth of the world’s popu- 
lation and the maintenance of basic subsistence levels for individuals 
mthin this population depend upon an accuratfe inventory of the world’s 
crops. Current proposed solutions to the problem of crop inventory 
depend on data gathered by satellites-remote sensing. Using ground 
reflectance patterns at selected wavelengths as data a discrimination is 
to be made among the various crops. Typically, to discriminate between 
fields of cotton and wheat a mathematical model is enployed that relates 
the remotely sensed reflectance patterns to their causes i.e. wheat 
canopies or cotton canopies , Many mathematical models are proposed by 
various investigators [2], [3], [5] with more or less unknown discrimina- 
tion capabilities. These models are roughly divided into two categories - 
deterministic and statistical. The deterministic models attempt to develop 
a cause-effect relationship between the input parameters (sun angle, viewer 
angle, plant parameters, etc.) and the observed changes in the reflectance 
patterns, with the ability to ej^plain observed changes in reflectance 
patterns. This type of model is desirable since the cause-effect relation- 
ship allows one to explain the alteration in reflectance patterns caused 
by crop diseases, pests, drought, and to optimize crop yields. The major 
weakness of this type of model is the inability to mathematically describe 
the si4)tle relationships that exists between the causes and effects. In 
fact this is such a serious problem that many researchers rely on stat- 
istical models using the philosophy that a crop can be identified by its 
reflectance properties by observing a large number of test plots of the 
same crop and measuring their variabilities. This type of mathematical 
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model has shmm good resuilts in laboratoiy environments but fails to 
discriminate adequately with a moderate diurnal variation [1] . 

Ihe puipose of this stutfy is to examine a proposed family of 
deterministic mathematical models for vegetative canopy reflectance 
developed by G* H. Suits of the Environmental Research Institute of 
Michigan (ERIWO [2] . A solution for each of these models is given, 
along with a computer program for inplementing these solutions. Pro- 
perties of these solutions will be discussed, and the mathematical 
models will be coiipared with actual field data. Finally, suggestions 
will be given for inprovements of these models and future areas for 
investigation will be discussed. 

2. The Solution of the Suits* Model for Canopy Reflectance 

The ijiportance of recognizing plant canopies by using remote 
sensing techniques has continued to grow. Mathematical models have 
been developed to interpret data acquired by remote sensing devices. 

It is hoped that mathematical models can predict different reflectances 
for different crops. A study of the models should reveal the weaknesses 
and strengths of each model. Also, a stu(fy of the different models should 
expose the features which can best differentiate between different crops. 

One of the models developed for identification of plant canopies 
was devised by Qvynn H. Suits. Suits' model marks the first attenpt to 
account for directional reflectance as a function of view angle. It 
also atteiipts- to trace changes in reflectance to spectral and geometric 
changes within the plant canopy. In Suits' model the plant canopy is 
divided into N layer which are infinitely extended. The last layer is 
always bounded by the soil. Each layer can have several conponents 
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(leaves, stalks, flowers) which exhibit different optical and physical 
properties. The con^ionents are assumed to be randomly distributed and 
homogeneously mixed. The con^3onents in each layer are i(fealized as a 
conhination of a vertical and a horizontal panel. The vertical and 
horizontal panels act as, Lambertian surfaces; that is, they diffusely 
reflect and transmit the incident light. The area of the horizontal 
panel is obtained by prdjecting the area of the conponent on a hori- 
zontal plane. Similarly, the area of the vertical panel is obtained 
by projecting the area of the coirponent on tiro orthogonal vertical 

planes . 

The radiant flux that interacts with the canopy is divided into 
two types: specular flux which arrives directly from the sun and 

diffuse flux. The symbols Ejj^(s,i,x), Ej^(+d,i,x), and Ej^^(-d,i,x) 
represent the specular flux, the upward directed diffuse flux, and the 
downward directed diff^e flux in the i^^ layer and level x for a parti- 
cular wavelength (x) . In Suits’ model the level of a layer is measured 
downward from the top of the layer. Since doivnward in Suits’ model is 
in the negative direction, then both the level and the total depth of a 
layer are negative nimbers. To determine ExC+d,i,x), E^(-d,i,x), and 
E;^Cs,i,x), the differehti^ equations 


dEj^C+d,i,x) 

dEj^(-d,i,x) 

dx 


- a. E,C+d,i,x)+b. E (-d,i,x) + c. E^(s,i,x), 
1 A 1 

a E,(-d,i,x)-b. E.(+d,i,x) - c’ E Cs,i,x), 

i A 1 


dE (s,i,x) 
A 

dx 


k. E^(s,i,x) 

X A 


and 
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must be solved. The constants a., b., c., c.’, and k are derived from 

1111 

measurements of the canopy conponents. If only one type of conponent 
occupies the i^ layer, then 

a. = Cl - t) a - S-iJL)], 

\ 

bi . = bh”hi= + Vv 

} 

=i' ' t°h”h ^ " 4’ Vv ® 1 ’ 

and 

■‘i = [°h\ * ® ^ 

where a is the average area of the projection of the component on a 
h 

horizontal plane, is the average area of the projection of the 
canopy component on two orthogonal vertical planes, n^ is the number of 
horizontal projections per unit volume, n^ is the number of vertical 
projections per unit volume, 9 is the polar angle for incident specular 
flux, p is the heinisperical reflectance of the component at this wave- 
length, and T is the hemispherical transmittance of the conponent at 
this wavelength. If there are more than one type of conponents in 
layer i, then the values a, b, c, c‘, and k are obtained for each type 
sei^arately and added together respectively to obtain a^, b^, c^, c,* , 

and k . For example, if there are two types of conponents in layer i, 
i 

\ 

then 

a^ = aCtype 1) + a(type 2) . 
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The boundary conditions require that at the top of the first 

layer the only do^vnward directed flux be specular flux. This can be 

stated as E. Cs,l,0) = 1. Hence the downward directed diffuse flux is 

zero, or E (-d,l,0) - 0. At the layer boundaries, the conditions require 
A 

that the uptvard and dcRmward directed flux be continuous . Finally , at 
the soil level, the boundary conditions require that all doimward 
directed flux be reflected to produce upward directed diffijse flux. This 
last condition may be stated as 

E^Cd,N,d^) = P3[E^(-d.N,d^) ♦ Ej(s,N,dj^)] 


where p is the soil reflectance, N is the last layer, and d^ is the 
depth of the last layer. 


2.1 Solution to the Boundary-Value Problem Associated with the N Layer 
Model ~ 


The system of differential equations may be stated as 

E.(x) = N. E- (x) for i = 1,2, ... , N, where 
i X 



■ dEj^C+<i,i,x) ' 
dx 


-a 

i 

b. 

1 

c. 

1 

i (X) = 

dE (-d,i,x) 
A 

s N. = 

-b. 

1 

a. 

1 

-c. ’ 
1 

dx 






dE (s,i,x3 
A 


0 

0 

k. 

1 


dx 




- 


■ Ej^(+d,i,x:> 
E^(x) = Ej^(-d,i,x) 

E^(s,i,x) 


For a given x, is the matrix representation 
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with respect to the natural basis of a linear transformation L of 
into r 2. If L has three distinct eigenvalues then the matrix representa- 
tion of L with respect to an eigenvector basis will be a diagonal matrix. 
The elements of the diagonal matrix will be the eigenvalues. In this 
case, the eigenvalues ,are,g, = (a.^ - b.^)V2, -g and k, . If 

^ ^ 1 X 1 X 

p + T = 1, then g^ = 0 and there are repeated eigenvalues. This is 
nearly the case in the infrared region which will be discussed in 
another section of this paper. On the other hand, if p + t < 1 then, 
with only one exception, there are three distinct eigenvalues. The 
exception occurs when 


9 = tan"i 


2 ) on ^ . m. 
" z. V v 1 

m- ^ 


where m^^ is the number of com- 


ponents in layer i. This angle of 9 causes = k^. 
An eigenvector associated with g^ is 


e, . = 
11 


a- + s- 
1 *■! 


1 


. An eigenvector associated with -g^ is 


®2i " 


a. - g. 
1 


1 


. Finally, an eigenvector associated witli 
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®3i = 


■li 


hi 


Ca. - k.) c. + c. ’ b. 
1 1 11 

a. 2 - k.2 - b-2 
111 


Ca. k.) c. ' + c. b. 
1 11 11 

a^2 - . 13^2 


Hence the matrix; 


representation o£ L with respect to S. = {el . , 

1 li’ 


21 - 


e^j^} is where 



g. 0 0 

1 

0 -g. 0 

0 0 k. 

1 


Now E. (x) = N.E. (x) 
1 11 


can be expressed as 


F Cx3 = D F.Cx) where F. Cx) = M.E. Cx) and F.(x) = M.E.Cx) where 
1 11 1 11 1 11 


M. = 
1 


a. + g. 
1 S] 


1 

S'. 

1 


a. 

1 


1 

^li 

- % 

1 

F. 

1 

^2i 

0 

1 


The solution to F . (x) = D.F. (x) is 
1 11 


FiM = 


A. e^i^ 


B. e'Si^ 
1 


C. e’'i== 
1 


where A^, and are constants of integration. 
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Thus 

E. fx) = M.F. (x) or 
11 

p.x g.x „ k.x 

E, (X) = A. eSi* e^. + Bj e^i e^. ♦ C. e x e^. . 

The constants of integration A., B. , and C must now be determined such 
that they satisfy the boundary conditions. 

At this point, a useful notation will be adopted to represent 
the boundary condition at the soil level. 


r 1 

a 


Ej^(+d,N,(^3 

b 

= a - p (b + c) . Thus f p 


c 


Ej^Cs 


will satisfy the last boundary condition. 

To determine A^, and C^, what will be done is to trade a 
boundary condition for an initial condition. Thus the boundary-value 
problem will reduce to an initial -value problem. 

Let Aq represent the upward directed diffuse flux at the top of 

the first layer. Then 



A 1 


r A, 1 


0 


1 

E^CO) = 

0 




1 


-Cl . 


det M. = C — * — 
1 a, + g, 


1 .) i 5 ^ 0, then M, is nonsingu3.ar. 


‘’i 
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Again since is nonsingular 
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■ % ' 


Ao 

= M 

B 

= M ... M, M 

0 

N 

N 

N N 11' 



Letting cl,, cl , and cl represent the coltaims of the 3x5 matrix 

... then = Ad boundary condition 

requires that £p(Ej^C<^)) = 0. 

Thus 

fpC%C%)3 = fpCAo clj clg) = Aq fp(elj) + fpCclg) - 0, or 


A = - — i— L. , It can be shown tliat for the one layer case 
^ £p(cl^3 

fp(clj^3 5 ^ 0. An attenpt was made to show 0 < Aq < 1 for the one layer 
case, but the attempt proved unsuccessful. Beyond the one layer case, 
it comes difficult to show that fp(cl^) r 0. Thus if fpCcl^] = 0, 
then this problem has no solution. However, if fpCcl^) ^ 0, then the 
constants of integration can be obtained from (1) . 

2.2 Solution to the Infinite Case for Suits' Model 

Another case yet to be considered is one in which, due to dense 
foliage or due to the depth of the canopy, the soil reflectance is 
negligible. In this case the boundary conditions become 

E;^(-d,l ,03 = 0, 

Ej^(s,l»0) =1, and 

ExC+d,l,x - «) = 0. 
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From the previous section, the solution 
ferential equations is (x) = e®i^ e^^ + 


to the system of dif- 

g.x ^ k-x 
e^i e. + C, e 1 e_. 

2 1 0 


Again, let A^ be the upward directed flux at the top of the single 
layer. Thus the initial condition for the infinite one-layered case 
becomes 


A 


^ ^ c 1 


' \ 

0 

0 


t g, - gi h 

1 1 r 


\ 

, 1 . 




- S . 



0 0 1 . 




Hence 


1 1 

■^1 b b ^1 ^2 ~ ^ = 1. The boundary condition requires 

11// 






that fp(Ej^Cdj)!) = 0. Thus A^ e ^ ifp(e^) + e ^fp(e^) + 


e^i‘^ifpCeg3 = 0, or fpCE^Cd^) = A^ e^®i‘^ifp(e^) -f- fp(ep + 

C S ^ d. 

e ^ ^ ifp(e^) =0. The depth of the single layer (d^) becomes 

large. Hence 


limit [A^ e^'i'^ifpCep + B^ fpCe^ + e^^i’^^i^^^ifpCeg)] = B^ fpCe^) = 0 

d^->-“ 

For 0 < p < 1, fpCe^^ ^ ^nd - 1, then 

ustag Aj i + B, Cj 5, = 0, 

(&! + hj) hi Cj’ + bj2 Cj 

A - . ■ — - — ^ ", • ■ ■ ■ 

^ b^+k^-a^ 

1 1 1 



12 


2.3 Solution to the Combination Case for Suits* Ivtodel 

The combination case was intended to describe the reflectance of 
a canopy primarily as initial surface reflectance from the top layer of 
the canopy. It is thought that the total reflectance from some crops 
such as ■\dieat originate mainly from light interaction with the top of 
the canopy crop (the heads of the wheat) . The combination case consists 
of two layers, the first having a finite depth and the second having an 
infinite depth. The flux densities of the first layer are described by 


EjCx) = Aj e®i^ e^i + e"®i^ e^j + e^^^ eg^ where A^, Bj, and 

are to be determined. Similarly, the second layer of the combination 
case is described by 


EjW - A, ej. 


+ B„ e 


•g2^ 


22 


+ C. 


k„x 


' 32 * 


The initial condition for the top layer is 


■ Ao ■ 


1 

1 



■ ^ ' 


ai + gi 

ai - Si 



0 


1 

1 




1 




^21 


- Cl . 



0 

0 

1 



Also A, - + B, - + C, £191= 0 which implies that B^ = -CAj + §21 

> bj ‘ b, 

In the second layer, the depth (d^) becomes large thus 

limit A2 e^^^^^fp (612^ ®2 + C2 e ^ ^ ^fpCe,^) ~ \ ^P^®22^ 

d2->“« 

This iiiT)lies that B2 = 0 since fp(e22) M for 0 < p < 1 . The boundary 

conditions require that E^( 0 ) = E^(d^), or 
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M. 




= M, 


B, 


'1 -J 


where M„ = 


^2 Sz 
1 


^2 ~ 

1 

b„ 


*12 


22 


and 




eel'll 

1 


el'll- 

a, gi 

C ^ 

^■1 " 



1 ■ 

-gl^l 1 

e’'A 


0 

0 



’ll 


’21 


Again, M, is nonsingular since det 1 


1 0 ^ 0- 


Hence 


■ Aa • 


> 

|W 

■ 

h 

= M -1 M, 
2 1 

B 

1 

= U-'^ M 
2 1 

- S - 


L ^1 - 



a2 + g2 ^2 - g2 l>2 

-CAi + 

1 


. C2) 


Let the elanents o£ the second row o£ be h, i, and j. Then 

Ajh - (Aj + %2i bj) i + 3 = 0» and 

£ ^5 i "1 * 

A ^ . I£ h = i, then there is no solution to this problem. 

1 h - i 

If h 5 ^ i, then A^, B^^, and Cj^ are known and A^, B^, and are obtained 
from (2) , 

Once the constants of integration have been determined for each of 
the three cases, the radiance due to an infinitesimal strip in the i 
layer can be calculated from Suits’ equation 
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4L, = [u. + V. ♦ „. SxC=.t.^J ] ^ 

liT iTr-* 

where u^ = I a^n^ i + (i] tan 0 

nij^ 

’"i ' ^ W’t Vv^ ^ 

“M. 

Wi = I O^n^ p + (i)2tan 0 tan 0 

"4. 

where 0 ib the polar viewer angle and m^ is the number o£ com- 
ponents in layer i. I£ i > 1, the £raction o£ radiance £rom this 
in£initesmal strip o£ layer i that is seen from the outside is 

i-1 - 

expC I % dj) expCf^x) where £^ = I aj^njj, + a^ny tan 0 and d. i 
j=l mj J 

the depth of layer j . The radiance as seen from outside the plant 

canopy becomes 

i-1 

ALj^ (outside) = exp( ^ £^ d.) expCfj^x) AL, . If i = 1, the fraction 

j=i J i 

sin^Jly becomes exp(f^x) and 


is 


aLj^ (outside) = exp(fjx) ALj^. 

The radiance of the entire canopy as seen from the outside becomes the 
integral of all such contributions plus the contributions of the soil. 

As an illustration, take the radiance of the infinite case. 

Now 

1 


Ej^C'*'<i,l,x) = 


a, + 


gi 
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E;^(-d,l,x) = ^ Aj and 


Ej^(s,l»x) = Thus 


ro 




£nX Ej^(+d,l,x) Ej, (-d,l,x) E, (s,l,x) , 

e ^ [u, + Vi -ii + Wi J dx 

'• i IT ■*■ 'll IT 


.^1 


Vl 


fo 


/-o 


[-4^ A * e’'l==] dx 


a + e, 
1 ®i 


A .. o8l^ 


knX. 


If. A: ® ' I *< 


ro 


+ 1 


£ix kix 

e e dx 


= "i [A. 1 


ai + gi £i + gi 


1 + C 1 — 1 

^ J? ^ J 


■1 f. ^k. 


V r. •*• J. _ -L 

* ji- 1^1 b. fl + gl ^2 fj * k. 


tr4Tr>- 

IT *1 


The refLectance (%) from the plant canopy then becomes 
h 

ik - „ FA £ £ + c £ 1 

E “ 1 + gj fl + ^1 

TA 1 1 , - 1 1 


+ w, [ 


1 > £^ + kj 


3. 
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Note that in the infinite case and combination case there is no contri- 
bution to the reflectance by the soil. In the N layer case, the con- 
tribution to the reflectance by the soil can be expressed by 

RCsoil) = e:qp[fj + f^ d^ + . . . + dj^3 E^C+d,N,dj^3 . 

2.4 Solution to the N layer Case for Suits* Model with Diffuse Light 
In the three cases discussed above, it was assumed that the only 
downward directed flux at the top of the first layer was specular flux. 
In cloudy days, however, a significant part of the light falling on the 
top of the first layer is downward directed diffiase flux. Also, the 
amount of downward directed diffuse flux at the top of the first layer 
in the visible region of the spectrum is a function of the time of day 
Csee Figure 1) as was shown by Jones and Condit [11]. Letting be the 
amount of do^vnward diffuse flux and be the amount of specular flux 
at the top of the first layer, then the initial conditions become 
E^(+d,l,o) = Aq 
E,^C- d,l,o) = 

E^(s,l,o) = Cq 
where ^ ^o ” 

The following discussion shows how the N layer case is affected 
by changing the initial conditions from 

^0 
Co 

A change in initial conditions will only change the constants of inte- 
gration A. , B. , and C. . Previously, the strategy used to solve for the 



wheie B and C„ are known values, 
o o 


SOLAR INCLINATION <DEG 
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FRACTIONAL CONTRIBUTION 

.0 *2 -4 < *6 1-0 


8H 




o 1 


o 

O 



CB 

O 


t 
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constants A. , B. , and C. was to solve for A first. It was found that 


1 1 


. fpCClg) 

A„ = ^ V . With A calculated, the constants A., B., and C, 

^ fp(Cli) o 1 -x’ 1 

were obtained from 


'A. 



1 


o 

B. 

= M.-1 M. , M. , ... M, M 

0 

1 

X X-1 x-l 1 1 




1 

1 




With a change in the initial conditions, the value of A^ changes. The 
value of Aq now becomes 

Bq fpccig) - Cq fpCCig) 

A = - . Again mth A„ calculated, the constants 

o fpCCli) ^ o 


B^, and are now obtained from 


F A. 1 


r A„ 1 

X 


0 

B. 

X 

•= M j ... 

Bo 

C- 





0 


The following changes in NTH2 (a computer program designed to yield the 
radiance from Suits* N layer model) were made to take into account 
changes in the downward directed flux at the top of the first layer. 

Lines 

1304 Print "Enter fractional contribution of both the specular flux** 

130fe Print **and Diffuse flux to the total downward directed flux."; 

1306 Input L3, L2 
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were inserted. Lines 1920 - 2000 were changed to 
1920 A2 = -1 * L2 A (G(l,2) -1 a SI * CGC2,2) + GC3,2))] 

1925 A3 = -1 * L3 * CG(1,33 -1 * Si a (GC2,3) + G(3,3))) 

1930 A4 = GC1,1) -1 A SI A (GC2,1) + G(3,l)) 

1940 IF A4 = 0 TOf 2480 

1950 A0 = (A2 + A3)/A4 

1960 R5(l) =0 

1970 R6(0) = 0 

1980 1(1,1) = A0 

1990 1(2,1) = L2 

2000 1(3,1) = L3. 

Similar changes were made on Inf 2 and Com 2 but are not included here 
since Inf 2 and Com 2 are special cases of NTH 2 (Inf 2 is the one layer 
NTH 2 with a large depth) . Tne resulting programs with these changes 
were called NTH 3, Inf 3, and Com 3. 

2.5 The Repeated Roots Case for Suits * Models 

For incident light in the wave lengths from 760 to 1250 n.m. 
the single leaf absoiption for a plant is small [4] . (See Figure 2) 

It can be shown that for small absorption, the Suits' differential 
equation coefficients a^^ and b^ are almost equal. If one assumes that 
^i ~ ^i eigenvalues of the matrix N^ as discussed in 2.1 are 
repeated, thereby yielding a different algebraic form for the solution 
to Suits' Model. Since this solution is relevant to the work done by 
Allen, Gale, and Richardson, (AGR ^fcldel) it is included to show the 
relationship betoeen their vegetative reflectance canopy model and the 
Suits' Model . 
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The matrix = 


-a. 

b. 

c. 

1 

1 

1 

-b. 

a. 

-c. ' 

1 

1 

1 

0 

0 

k. 


has characteristic 


polynomial c(X) = giving eigenvalues 0, 0, and k^. An eigen- 

vector corresponding to the eigenvalue k^ is 

^ I 

I ^ 1 H t • 

( 1 ) 


e = 
3 


Where = c-(a- - k^) - c^» 


p _ b. c. + c. ' (a. + k.) 
_ 11 1^1 l'^ 

k.2 

1 


( 2 ) 


To find eigenvectors (if they exist) that form a basis for the 
subspace N? x = 0 of three -space we must find non- zero vectors e^^ and e^ 


such that 

e^ = 0, ®2 ~ ®i* 




'1 


U «3 J 


then 


ai Oj ♦ b. ♦ Cl “3 = 0 


b. a + a. a - C. a = 0 
i 1 12 13 


k 0 , = 0 


a. must be zero, and a solution is a, = 1, a„ = 1 so that 




1 

1 


0 
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If a,- 



then the condition that 


N e = e, forces 
1 2 1 


- * ‘=1 ®2 ^ ‘ ®i 

- bi 6i + a. Sj - Cj'Bj = a. 


k. B, = 0 
1 3 


A solution is = 0, 8^ “ 
'-1 

= 0 . 

0 


The set S = {e^, e^, e^} is a basis for three-space, 
representation of with respect to S as Since 

Mg =0, Ms e^ = e^, e3 = k. eg 



■ 0 

1 

0 


’ X 

1 

“s = 

0 

0 

0 

. If we let X = 

X 

2 


0 

0 



X 

3 


becomes 

^3 “ \ ^ 3 ’ 

yielding 

lc*x 

X = C e ^ X = B , X = B X •!• A, 
Xg - e -^-2 1’ 1 1 1 

where Ajj B^, ars constants of integration. 


Denote the 


the DE system 




4 
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The general solution to the original system is thus 


B = Xi e, t Xj * Xj 


k.x 


= CAi + 3^x) e^ + Cg + e 1 e 


so that 


Ic-x 


E(+d,i,x3 = (Aj - Bp + B^x + 5^ e i' 

E(-d,i,x) = A^ + fljX + Cj 5^ eV 
k.x> 

ECs,i,x) = Cj el.. 

For the one layer case id.th the initial conditions 


A 

c 

0 

1 


it is seen 


that = 1> Aj = - ■ ^2 " = 1 “ P. 

£p fe ■) = -1 , fp fe ) = ? - p CC +1) then the boundary condition gives 
2 3 1 2 

at X = -d gives 0 = - Ap(-d)(l-p) + ~ 5^ " V M) 

+ e'^i*^ - pCC^ + 1)], or 

A = (?! - S;)dCl - P) Cl - e'V [gi - P(S2 * m 

° 1 + d(l - p) 

yielding a solution to the one layer finite case. For the one layer 
infinite case, let d-j--<», which gives A^ so that = 1, 

= 0, A, = - 5^. ■ 

The infinite solution is 


I 


f 
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E(+d,l,x) = -5 +5 
2 1 

E(-d,l,x) = eV 

ECs,l,x) = 

The N-layer repeated problem is solved in a manner analogous 
to the N-layer non-repeated case discussed in 2 . 1 . 

1 -1 Cj ' 

IOC 

2 

0 0 1, 

with C • and C - the eigenvector con^onents using the parameters of 
X 1 2 ^ 

layer i in equations (i) and ( 2 ) . The unknown initial condition 
solves the equation 

A^fpCGl^) + fpCClg) = 0 where [Cl^, Cl^, CI3] = 

M = 

columns of M respectively) . 

If fp(Cl^) 0 , the system has a solution with 

A - - fpCClg) 

° ~ fpCtli) ' 

To find a relationship betvreen the Suits and AGR Models one 
uses the one layer Suits Model \ifith repeated roots. Let = b^, 
(which is equivalent to asstmiing that y = 0 for the AGR parameter for 
single leaf absorption). The following table i^elates the Suits para- 
meters to the AGR parameters . 




' 1 

-Ci+d.) 

^li 

-k.d. 
e 1 1 

M. = 

1 

-d. 


-k.d. 
e X 1 

1 

0 

1 

0 

21 

e-^i'li 
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Suits Parameter 

AGR Parameter 


y * + 3 ' + p* 


IJ + 3 

b. 

3 

1 


S’ 

F' 

c. 

3' 

1 

] 


0 

? 

EC+d,i,x) 

sCn) 

E(“d,i,x) 

t(n) 

E(s,i,x) 

Un) 

Suits Conditions 

AGR Conditions 

E(-d,l,0) = 0 

t(0) = 1 

E(s,l,0) = 1 

1(03 = 1 

E(+d,l,r) = p[EC-d,l,r)+ ECs,l,r)] 

s(r3 = 0 

positive direction upward 

positive direction dowrarard 


Table 1. Analogy Table for Suits and AGR Models 

Using this table, it can be seen that the differential equations 
describing the t\'/o models are identical , but the boundary conditions 
are different. 

As the AGR Model was not designed to measure light reflectance 
exterior to the canopy , the look angle 0 in the Suits Model does not 
have an analogy. Experimental evidence [5], [6], indicates that these 
models agree mth actual data taken within a plant canopy . 
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The AGR parameters were calculated by use of regression analysis 
on e^qjerimental data. With the AGR coefficients calculated for a corn- 
field in Ithica, New York [5] an attempt was made to calculate the Suits 
coefficients, but the conditions inposed on the Suits coefficients by 
Table 1, forced several of the Suits Parameters ' to be negative, a con- 
dition that is physically impossible. ' 

3. Theoretical Light Penetration in a Plant Canopy Using the Suits 
Model 

This section of the paper uses the Suits Model to discuss the 
inportance of canopy depth in reflectance measurements. These deter- 
minations were based on the following definition: 

the penetration depth of a plant canopy at a given wave- 
length of light is that depth at whici reflectance 
obtained from the one layer finite Suits Model is within 
5% of the reflectance obtained from the infinite Suits 
Model. 

Figure 3Ca3 illustrates a determination of the penetration depth 
for eoperimental data collected for a cotton canopy at 850 n.m. As 
the canopy depth increases, the reflectance increases due to the in- 
creased number of leaves, each leaf acting as a good reflector at this 
wave length (the absoiption of a single cotton leaf is .018 at this 
wave length) . The reflectance is within 5% of the infinite canopy 
reflectance Cthe horizontal asynptote) at a depth of 47 cm. , indicating 
that vegetative structure below 47 cm. in the canopy contribute less 
than 5% to the total reflectance at 850 n. ra. Thus, only 47 cm. of the 
cotton canopy is sanpled by 850 n. m. light. That is, 850 n. m. light 
can be measured as a conponent of the total light transmitted to lower 
depths in the canopy, but would fail to exit the canopy tpward into the 




DEPTH <CM 


27 


REFLECTANCE C^L/E) 



28 


field of view of the observer. The behavior of 850 n.m. light is 
characteristic of most infrared radiation with respect to penetra- 
tion depth. 

Figure 3(b) indicates an interesting contrast that occurs for 
650 n. m. li^t. As the canopy depth increases, the reflectance 
decreases, due to the increased nimiber of leaves. Each leaf is now 
acting as a good absorber at this wave length (the absorption of a 
single cotton leaf is .877) . Reflectance decreases to within 5% of the 
infinite canopy reflectance (the horizontal asymptote) at a canopy 
depth of 16.4 cm. This behavior is characteristic of roast light in 

the visible spectrum. 

The piocess of calculating the penetration depth for various 
wave lengths was continued, and the results are shown in Figure 4. 

To illustrate the difference between the visible and infrared pene- 
tration depth, the convention was adopted of assigning a negative 
sign to the penetration depth if the reflectance decreased to its 
limiting value and a positive sign if the reflectance increased to 
its limiting value. This graph indicates a greatest depth of 47 cm., 
so that only the first 47 cm. of the plant canopy are useful for 
reflection measurements at any wave length. As expected, infrared 
radiation penetrates to the greatest depths in the cotton canopy and 
the visible wave lengths penetrate the canopy to the least depth. 

Of great interest is the penetration depth for an estimated 695 n. m. 
radiation (no data was available at the wave length, and this result 
is based on interpolation) . At this wave length the graph indicates 
'that observed reflectance data fails to penetrate the canopy, and what 
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Figure 4 

The Penetration Depth of Light in a Cotton Canopy as 

a Function of the Wavelength. 
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is recorded corresponds only to surface reflections. In sumnary. ^or 
the 189 cm. cotton canopy considered, incident light penetrated to a 
depth no greater than 47 cm. indicating that the soil reflect^ce data 
was not utilized in the model. 

Allen and Richardson [7] have observed phenomena analogous to 
this with single leaf models. For the cotton crop considered = 104, 
nj^ = .00125 and “ *13. For a maximum depth of 47 cm., tlie 
cuirailative L.A.I. is .13 x 47 = 6.11, so that infrared light penetrating 

4 

to this depth and exiting the canopy to the observer must pass through 
12 leaves, and by an analogous calculation, light in the visible region 
must pass through 4 leaves, Allen and Richardson found that by stacking 
cotton leaves in a spectrophotometer, the reflectance ceased to vary 
mth 8 cotton leaves in the infrared region and 2 cotton leaves in the 
visible spectrum, giving approximate agreement with the results found 
by use of the Suits ^fodel. 

On the basis of this limited data, one can conclude the following 
rules for use of the Suits Models: 


(a) Since the maximum penetration occurs in the infrared, the 
effective depth d of a canopy is no larger than d = b/aj^n^^, 

(b) If the canopy to be considered has depth greater than d, 
use the infinite model, and 

(c^ if the canopy considered has depth less than d, use the 
finite Suits Model and collect ground reflectance data. 

These relationships are summarized in figure 5 for a typical plant from 


a canopy. 

Penetration depth is a useful concept in explaining canopy 
reflectance for crops with layered structures such as wheat. As an 
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Figure 5 

Relationship between Penefation Depth and Wavelength for a Typical 
Canopy. 

a. 600-690 n.m. with small penetration into the first layer. The 
1 layer infinite Suits' Model fits the data well. 

b. 500-600 n.m. with moderate penetration into the second layer. 
The combination case for Suits' Model fits the data well. 

c. 690 n.m. and above with deepest penetration. A multi-layered 
finite Suits' Model fits the data well. 
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exanqJle, consider Figure 6. In this case, two different Suits Models 
were employed on wheat data, the first used an infinite Suits Model , 
and the second a four layer finite Suits Model. The first model used 
only the parameters collected for the wheat heads and the second model 
used all parameters collected for each layer and the soil reflectance. 

As would be expected, the too graphs agree in the visible region, as the 
penetration depth is small and the second model is utilizing only the 
data from the wheat heads. The difference betoeen the too graphs become 
noticeable around 700 n. m.r, where light has penetrated the layer of 
wheat heads in the second model and begins to sample the dissimilar 
reflectance properties of the green leaves. This difference continues 
to change until around 800 n. m. , where the penetration depth has 
reached the soil level, and remains essentially constant thereafter. 

This exanple stresses the fact that reflectance data from around 650 n. ra. 
to 850 n. m. contains a great deal of information on the internal 
reflective structure of a plant canopy. 

4. The Interchange Property for Suits Models 

This portion of the paper will establish a rather non- intuitive 
property of the Suits Nfodel - the interchange property. In general, 
for an isotropic plant canopy, the reflectance is a function of the 
sun angle 9 and the observer angle 0 measured from the vertical, and 
can be written as 

canopy reflectance = R(9,0]. 

The interchange property thus states that 


RCe,0) = 
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or that the same reflectance results whenever the position of tlie 
observer and the sun are interchanged. 

The property will be established by use of tlie infinite case 
of Suits Model for non-diffuse light. From section 2.2 of the paper, 
using the infinite model, 


RC9,0) = {■ 


-b[bc - 1 - Ca h)c'] 

(g2 - lc2)(a + g)Cg + k'3 




c(a - h) + bc» 

(g2 - + kO 


u 


+ (a + k)c *3 

(g2 . kZ) (g + k’ ) 


be Ca + k)c^ j y 
CgZ - kZKk + k‘) 


(13 


+ (: 


k + k> 


} w. 


where these parameters are ^defined in section 2.1. One can observe a 
good deal of symmetry in tliese parameters, in fact, if one interchanges 
the roles of 0 and 0 , k becomes k' , k' becomes k, u becomes c* , c' 
becemes u, v becomes c, c becomes v, and w remains unchanged. 

Expression (1) is formidable and the key step in proving this 
result is to re-write this expression. After considerable trial and 
error, one re-wites ( 1 ) as 

R( 9 , 0 ) = -b[bc + (a - g)c* ] u - [be (a - g)c*]v (a g) ^ 

2gCa + g3 Cg + k) Cg + k*3 

{ [bc+Ca+g3c» ] [bu-fCa-j-gpvIKg-i-kj-i-Ca-^g) (gfk'}{ [cCa+g3-<-bc 3m[bc->-(a-g3c» jv) 

2 g(a+g) Ck+k' 3 (g+k * 0 Cg+k) 


+ {■ 


k + k‘ 


w. 
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V/hile the first two terms are quite complicated, it can be seen 
by observation that both denominators are symmetric in k and k’ , thus 
one needs only show that the two numerators are symmetric in 9 and 0. 

I ' 

The third teim is already symmetric in 0 and, 0 , , 

Define the first numerator as 
Nj(0,0) = -b[bc + (d-g)c']u - [be + (a-g)c*]V Ca+g) 
then * 

N^(0,0) = -b[bv + Ca-g)u]c' - [bv + (a-g)u]c Ca+g3 
or, rearranging 

Nj(0,0) = [-b(a-g)c’ - cCa+g)(a-g)]u + [-b^c’ - bc(a+g)]v. 

Using the fact that b^ = •* g^ and factoring gives 

NjC0>0) = -b[bc + Ca-g)c']u - [be + Ca~g)c']v (a + g) = NjCQ,0). 
In a similar manner, define the second numerator as 
N 2 CS> 0 !) = {[be + Ca+g)c'][bu + (a+g)v]}(g+k3 + 

(a+g3 Cg+k*){[c(a+g3 + bc’]u + [be + (a-g)c']v}. 

N2(0»9) = {[bv + (a+g)u][bc‘+ Ca+g)c] ICg+k’) + 

Ca+g3 (g+k{ [vCa+g3 + bu]c' + [bv + Ca-g)u]c}. 

These terms can be. re-arranged so that 

N2(0,9) = {(a+g)u[bc' + Ca+g)c] + bv[bc’ + Ca+g)c]} (g+k') + 


(g+k){[v(a+g)2 + b(a+g)u]c’ + [b(a+g3v + Ca^-g2>u]c}. 
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But = a2 - and factoring gives 

N2(0,9) = {u[bc* + Ca+g)c] + [(a-g)c' + bc]v}(a+gKg+k’) + 

(g+10{[bc + Ca+g)c’][bu + (a+g)v]} = 

Thus, it has been shown that 
R(9,J3) = RC0,0). 

Considering that the Suits Model is a generalization of a math 
model used to describe light reflection from a stack of glass plates, 
and that in most optical models the interchange is valid, then this 
result should not appear too surprising. 

It should be noted that the interchange is not valid in the 
presence of diffuse light i.e. E_<j(0) > 0, EgCO) > 0, E_d.(0) + Eg(0) = 1. 
(See 2.4) 

For example, in the case that E-dCO) = -6, EgCO) = .4 
R(10O,50O) = .752371, but R(50O,10O) = .736506. (The infinite diffuse 
light Suits Model applied to the parameters for cotton in data set 1.) 

Both the N- layered Suits Model and the Combined Suits Model 
exhibited the interchange property for non-diffuse light for all data run 
on the computer, but due to the complexity of the equations involved for 
reflectance, an analytic proof of the interchange property was not 
attenpted. Similarly, both of these models did not exhibit the inter- 
change property for diffuse light. 

5 . Directional Reflectance for the Suits Models 

The general problem of calculating the canopy directional 
reflectance as a function of the sun angle, viewer angle, and viewer 
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azimuth has been reduced by the Suits assunption o£ isotropic canopy 
conditions mth respect to azimuth, so that the reflectance is a 
function of only the sun angle and viewer angle measured from the 
vertical. It is ej^ected that such a simpli:fying assunption might lead 
to errors between experimental data and theoretical calculations of the 


bi-directional reflectance function. 

This section of the paper will discuss the general shape of the 
surface R(9,^) (the bi-directional reflectance function) and derive 
several useful foimulas. Using the N layer Suits* Model, the radiance 
contribution from the i"th layer is given, from Suits equations as 


R. 

1 


P. 

1 


rxi- 


X. 

1 


h.X .r' J 

e 1 ALj^dx 


where the i"”*^ layer has boundaries x^ and x^_^ respectively, 

i-1 

r \ 

'■ i 

y 


X X 

. = 1 if i = 1, and P = exp( I if i > 1. 

1 1. . i=i J J 


P. = 
1 

But 


AL^ = uE+^ + vE.jj + wEg where 


“ “ ”h"h ^ * Vv 
■ ■’h”h ♦ Vv ® 


w 




and E^^ = . e‘V + . e^i^ + C^. e^i^ 

E . = e‘Si^ . B„. e-i^ . C„. 
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21 


E = C eV 

Si ® 
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Observe that A. B.., C.. for i = 1,2,3, j = 1,2,3 are independent 
ij ij 

of the view angle 0. For i = 1, (the first layer) = 0. Inte- 
gration and re-arrangement gives 

_ uAjn + VA21 + VE21 uCji + VC21 + WC31 

R, + — + — — 

^ k - » k + I, k + k. 


- ^^21 - expCCk, - g,)x ) 

K, - g, 


^11 expCCk. + gjx.) 

Ki + gi 


uC,. vCg^ •>• wCg^ 

K + k 
1 1 


expCCK^ + 


As 0 ^ I. ^ each of tlie last three terms involving the exponential tend 

P T 2 

to zero, and the first three terms approach — ^ — [E(+d,l,0) + tan 9]. 

For i > 1, as 0 ^ , each of the seven terms resulting from 

integration approaches zero. Thus -»■ 0 as 0 y for i > 1. 

The ground radiance contribution after N layers is given by 
Pj^(+d,N,x^), and since 0 as 0 j and E(+d,N,Xj^) is' 
independent of 0, then the- ground radiance contribution goes to zero. 

Thus, we have shown, that as 0 approaches j , the only contri- 
bution to the reflectance is from the top surface of the canopy . In fact , 

\ 

Lim RC0,0)= ^-|p [E+a<^0) + I tan 9] . CD 



This result indicates for 9 and 0 sufficiently close to , 
R(0,0) becomes larger than one, a result that is physically nnpossible. 
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For non-di££use light a typical surface for a single wave length 
of light, as calculated by the conputer forms the shape of a saddle. 

Figure 7 is an exan^ile of such a surface for a cotton canopy with the 
iso-reflectance lines sho^vn to represent the surface contours. One can 
observe that the lapper right comer of the graph becomes larger than 
one as indicated by (1) , and that the remainder of the graph exhibits 
only a small amount of variation. 

For diffuse light of a single wave length with a proportion of 

\ 1 

m direct light and 60% diffuse light the reflectance surface for a 

cotton canopy is shown in Figure 8. 

In this case one can observe that the interchange property fails 
to hold, and the contours of Figure 7 erode into the contours of Figure 8. 

6. A CompaTison of Experim ental and Theoretical Results for Plant 
Canopies 

A con?)arison of Suits’ Reflectance (^Ve) values and the 
experimental values for cotton at various wave lengths for the same 
sun and viewer angles (Figures 9a and 9b) shows that the qualitative 
trends exist. The maxima and minima do not occur at the same wave 
lengths due to a calibration problem in measuring single leaf reflect- 
ance values with the Beckman DK- 2A spectrophotometer . The calibration 
problem caused a shift in wave length of about 30 nm. As can be noted, 
the theoretical values are significantly greater than tlie experimental 
values in the infrared region. Quantitative comparisons between the 
theoretical and experimental values are not possible at this time for 
several reasons. Among those reasons are normalization problems 
betiifeen the theoretical and experimental values. Also, there are varying 
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Figure 7 

Iso-Reflectance Lines for the Reflectance Surface 
R(0,0) of a Cotton Crop. 0 is the Sun Angle and 0 
is the Viewer Angle. 
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Figure 8 

Iso-Reflectance Lines for the Reflectance Surface 
R(0,0) of a cotton crop using 40% Specular Light and 
60% Diffuse Lights. 
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Figure 9 

Theoretical and Experimental Reflectance Curves for 
Wheat and Cotton as a function of wavelength. 
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methods used by ejqperimentalists to describe directional reflectance 
not all of which yield equivalent results . Another problem is the 
finite angular field of view that experimentalists must cope with to 
collect enough light for meter readings. In theoretical models, an 
aspect in question is the factor -rr that converts radiance to flux 
density. This is a well-established relationship for lambertian sur- 
faces. It is not clear however that n converts from radiance to flux 
density for non- lambertian surfaces such as canopy surfaces. 

Results similar to those of cotton were obtained in conparing 
Suits* Reflectance values and the experimental values for wheat 

(Figure 9c and 9d) . Again there is a shift of about 30 nm. between the 
experimental and theoretical results due to calibration problems. 
Similarly the theoretical values are greater than the experimental 
values in the infrared region. A conparison of Figures 9a and 9c shows 
that Suits* Reflectance (^Ve) values fo-^ ton are significantly 
greater than the Reflectance values for wheat. This is 

primarily so because of the different absorption properties of the 
wheat compc'-.ents in the infrared region and also because of the soil 
reflectance values which play a significant role in canopy reflectance 
for wheat due to the low LAI of wheat. 

7. Conclusions 

Both the Suits and the AGR models for vegetative canopy reflect- 
ance are generalizations of the Kubelka-Munk theory for the trans- 
mission of dif fiise light [8] . Both models assume incident specular 
light but assume that this light, upon being transmitted through a 
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leaf is diffijse enough to satisfy the KubeHca'-Munk hypothesis. In the 
words of [5] "Although a typical leaf such as com has strong specular 
reflection, the structure of the leaf is such that the transmitted 
light is very diffuse". This article cites [9] as an example of the 
strong directional reflectance properties of leaves. Both of these 
models, as experiments veriify, model light intensity interior to the 
canopy well, but fail to accurately model light reflected exterior to 
the canopy. [10] indicates experimental directional reflectance trends 
for a wheat field not found in the Suits Model. This deficiency in 
both Kubelka-^^unlc Models is due, in the opinion of the authors, i to the 

I 

failure in consideration of the strong specular reflectance from the top 
surface of the canopy. An initial atteirpt was made to alter the Suits 
Model by the development of the combined case, using a thin reflective 
top layer, and an infinite second layer, vdiich resulted only in a 
partial success. On the top surface of the canopy, one must consider 
both the plant geometry, distributions, axial "droop", etc.) and 

the single leaf directional reflectance as measured by [9] , rather than 
the single leaf hemispherical reflectance. The transmitted light 
through the top leaf surface is diffuse and a Kubelka-Muhk type model 
such as the one proposed by Suits should be adequate to describe light 
reflected from mthin the canopy. Thus, it is felt! that surface specular 
reflectance is a phenomena that must be modeled, and can be adjoined to 
the existing Suits Model as a refinement. 

The existing Suits Models, on the other hand, showed qualitative 
agreement behveen the experimental results gathered for cotton and wheat 
and the theoretical calculated results for reflectance versus wave length. 
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Again, the authors felt that qualitative agreement could be converted 
to quantitative agreement either with a revised analysis of the con- 
version factors necessary to convert bettveen radiant intensity and 
radiance for a non-Lairibertian source or a re-evaluation of the experi- 
mental proceedures needed for normalization of field data. 

The Suits Models agreed with work done by Allen and Richardson 
on the infinite reflectance of stacked leaf models, indicating that 
light will penetrate only through canopies with a cumulative LAI of less 
than 6, and that infrared light penetrates vegetation better than light 
in the visible region. The Suits Model gives a quantitative way of 
discussing light penetration through a canopy by means of the penetra- 
tion depth. The model indicates that light incident on a plant canopy 
is reflected to an external observer from very near the canopy surface 
in the visible region and the reflections progress to vegetative deeper 
within the canopy as the light wave length progresses into the infrared 
reaching a maximum depth at about 850 nm. and leveling off thereafter. 
Thus, the authors believe that in the wave length region from 650 n.ra. 
to 850 n.m. light scans a vertical profile of the plant canopy , and 
most of the information about the plant canopy is contained in this wave 
length region. This hypothesis is especially useful in explaining the 
differences between layered crops such as wheat and homogeneous crops 
such as cotton. 

It is felt that continued work on the Suits Models should pro- 
gress in the areas of data normalization, for both the experimental and 
the theoretical. Further, a "hybrid" model should be constructed that 
incorporates the specular directional reflectance of surface leaves and 
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uses the Kubelka-Munk theory for the diffuse light that exists below 

the top surface. Both a lack of time and a lack of experimental 

data on both single leaf and canopy directional reflectance prevented 

the authors from completing such a model. Finally, the authors feel 

that work should continue on a study of the differential coefficients 

proposed by Suits. Some of the coefficients could be derived by the 

authors, but others were in direct conflict both with the author's 

results and established facts. But, since the purpose of this paper 

was to test the Suits ^fcldel, all coefficients used were those supplied 

by Sui.ts [2]. This choice of coefficients explains, for example, why 

the reflectance is greater than one for large observer and sun angles 

on the directional surface R(9,0) pictured in Figure 7, and could 

possibly explain the relative insensitivity of the model to changes in 

the parameters a, n and a n . 

h h V V 
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Appendix 1 

A brief discussion 6n how to obtain the theoretical graphs of 
this paper follows. The experimental graphs were obtained with data 

I 

supplied by Dr. Edwin LeMaster. Since the total LAI(aj^nj^dj) of the 
cotton crop was 24.6, the one layer Infinite case was generally used 
on cotton. The four layer N, layer case was used on wheat. 

Figure 3. Graph a ^las obtained using Data set 1 (with 8 = 0°, 

0 = 0°, and X = 850 nm.) and using the Suits’ one layer model, using 
the con^Juter program NTH2. Graph b was obtained using Data set 1 
(with 0 = 0°, 0 = 0°, and X = 650 nm.) and using the Suits one layer 
model, using con 5 )uter program NTH2. The horizontal asymptotes are the 
Infinite case reflectance value at the Wo wave lengths respectively. 

Figure 4. The penetration depth at a particular \^ave length was 
obtained using Data set 1 (with 8 = 0°) and iising the conputer pro- 
gram NTH2. The depth was varied until the Reflectance (^Ve) using the 
conputer program NTH2 was within 5% of the Reflectance (^VE) using the 
conputer program Inf 2 for a given wave length. 

Figure 5. This graph was obtained using Data set 2 (mth 8 = 15° 
and 0 = 0^) and using Suits* four layer model, using the computer pro- 
gram NTH2. 

Figure 6. Gr^h a was obtained using all of Data set 2 (with 
0 = 15° and 0 = 0°) md using Suits' four layer model, using the com- 
puter program NTO2. Graph b was obtained from the data on heads in 
Data set 2 ( with 0 = 15° and 0 = 0°) and using Suits' one layer model, 
using the computer program Inf 2. 
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Figure 7. This gr^h was obtained using Data set 1 (with 
\ ~ 850 nm.) and using Suits* one layer model, using the conputer 
program In£2. 

Figure 8. This graph was obtained using Data set 1 (with 
X = 850 nm.) and with E (s,l,0) = .4 and E (-d,l,0D = .6 using Suits’ 

A A 

one layer model, using the computer program In£3 which is th§ altered 
In£2 program. 

: ! 1 •• 

Figure 9. For both the experimental and theoretical graphs in 
Figure 9, 0 = 15° and 0 = 0°. Graph a was obtained using Data set 1 
and using Suits* one layer model, using the computer program Inf2. 

Graph c was obtained using Data set 2 and using Suits* four layer model, 
using the computer program NTH2. 
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Data Set 1 
Cotton 


Layer 1 13 green leaves 

a, = 104 cm. 2 
h 

= 85 cm. 2 


X(nm.) 


500 

.099 

.107 

.100 

550 

.111 

.145 

.126 

600 

.076 

.078 

.160 

650 

.073 

.050 

.175 

700 

.308 

.382 

.190 

750 

.443 

.531 

.325 

800 

.444 

.540 

.417 

850 

.442 

.542 

.396 

900 

.438 

.542 

.398 

950 

.430 

.539 

.380 

1000 

.434 

.547 

.342 

1050 

.432 

.551 

.353 

1100 

.422 

.546 

.389 


nj^ = .00125 1/cm. 3 = 189 cm. 

n = .0007 l/cm.3 


Bata Set 2 
Wheat 


Lay 'r 

1 

1) heads 



depth = ,09 m 



= .000369 


6170 l/m^ 




= .001015 

n = 
V 

6170 l/m3 


Layer 

2 

1) green leaves 



depth = .20 m 



= .00144 , 

^h = 

4360 l/m3 




cT = .00107 

V 


4360 l/m3 




2) stems 






a, = 0 
h 


0 l/m3 




0^ = . 0006 

n = 

V 

2780 1/m^ 


Layer 

3 

1) green-broivn 'leaves 


depth = .10 m 



a. = .00148 
h 


8330 l/m3 




= .000872 

n = 

V 

8330 l/m^ 




2) steins 






a = 0 
h 

"h = 

0 l/m3 




= .0003 


5550 l/m^ 


Layer 

4 

1) broivn leaves 



depth == .12 m 



a, = .000259 
h 


12350 l/m3 




= .0000646 
2) steiiis 


12350 l/m^ 


, 


a. = 0 
h 

^h" 

0 l/m^ 




CT - .00048 

V 


43o0 1/m^ 



Green 

Efl Heads Green Broivn 

leaves ■ leaves 

X p T p X p X 


500 

.126 

.005 

.115 

.051 

.154 

.017 

550 

.137 

.004 

.119 

.066 

.177 

.021 

600 

.115 . 

.003 

.095 

.033 

.146 

.013 

650 

.110 

. 0-03 

.092 

.027 

.136 

.021 

700 

.371 

. 005 ' 

.308 

.258 

.366 

.169 

750 

.578 

.011 

.453 

.381 

.452 

.266 

800 

.591 

.023 

.453 

.389 

,457 

.306 


850 .595 .038 .452 .392 - .460 .333 

900 .577 .041 .446 .394 .456 .346 

950 .543 .036 .434 .394 .450 .355 

1000 .571 .053 .439 .400 .453 .362 

1050 .572 .061 .437 .402 .449 .367 

1100 .530 .047 .423 .398 .440 .368 




1 - 


Brown Stems 

leaves 

p TP Ps 

.203 .007 .115 .004 .129 

.243 .008 .124 .004 .160 

.287 .013 .093 .003 .176 

.329 .021 .086 .003 .222 

.370 .031 .367 .035 .203 

.402 .043 .595 .102 .294 

,433 ,057 .605 .113 .293 

.454 .069 .605 .114 .306 

.472 .077 .589 .103 .329 

.486 .086 .552 .080 .348 

.497 .094 .579 .102 .397 

.504 .101 .583 .109 .376 


509 


104 


546 


081 


333 


I 
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Appendix 2 

This section contains a listing of the computer programs 
associated with the three cases of Suits' Model discussed in this paper. 
The conputer program NTH2 was designed to calculate the radiance of 
the N layer model of Suits. Similarly, Inf 2 and Com2 were designed to 
calculate the radiance of the Infinite case and Combination case 
respectively. Note that the computer programs yield the radiance and 
not the reflectance from a plant canopy. To convert from radiance to 
Reflectance (^Ve) , the radiance must be multiplied by it. 
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NTH2 10;18 27-MAY-75 


1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 


PRINT'DD YOU WISH TO READ INSTRUCTIONS < YES OR NO) '5 
INPUT A$ . t • 

IF A$0*YES* THEN 1230 ^ 

PRINT' THIS PROGRAM DETERMINES THE RADIANCE* AS PROPOSED BY 
PRINT'GWYNN SUITS* FROM A PLANT CANOPY CONSISTING OF N LAYERS.' 
PRINT'THE PROGRAM WAS WRITTEN IN PARTIAL FULPILLMENT OF THE THESIS* 
PRINT'REOUIREHENTS Ol=" JUAN H, CANTU* WHOSE GRADUATE ADVISOR WAS* 
PRINT*DR. JOE CHANCE.* 

PRINT* THROUGHOUT THE FOLLOWING CALCULATIONS* THE SAME UNITS OF 
PRINT 'MEASUREMENTS SHOULD BE USED. THAT IS* IF THE DEPTH IS IN' 
PRINT'CM.r THEN THE AREA IS IN CM. **2* AND THE VOLUME IS IN CM.^r*3. 
PRINT* WITHIN A SINGLE LAYER THERE CAN EXIST SEVERAL COMPONENTS* 
PRINT* (LEAVES* STALKS* FRUITS* FLOWERS) EACH WITH ITS OWN SET OF* 
PRINT’PHYSICAL AND’ OPTICAL PROPERTIES. THE AREA OF A COMPONENT* 
PRINT'IS DIVIDED INTO THE AREA PROJECTED ON A HORIZONTAL PLANE AND* 
PRINT’THE AREA PROJECTED ON TWO ORTHOGONAL VERTICAL PLANES. IK* 
PRINT'THE AREA IS MEASURED IN CM.^*2* THEN THE NUMBER OF HORIZONTAL* 
PRINT 'AND VERTICAL COMPONENTS PER UNIT VOLUME SHOULD BE IN 1/CH.#*3.* 
PRINT'THE TRANSMITTANCE <T) AND REFLECTANCE <R) OF A COMPONENT ARE* 
PRINT 'UNITLESS AND 0<T J * 0<R<1» O-T+Rll. THE SUN ANGLE AND VIEU** 
PRINT'ANGLE ARE MEASURED IN DEGREES. BECAUSE OF THE CHOICE OF CD-* 
PRINT' ORDINATE SYSTEM* THE DEPTH OF A LAYER SHOULD BE ENTERED AS* 
PRINT'A NEGATIVE NUMBER.' 

DIM E(3*3)jI(3*1> 

PRINT'ENTER WAVELENGTH IN NANOMETERS'* 

PRINT'ENTER NUMBER OF LAYERS* BUN ANGLE* VIEW ANGLE* AND REFLECT-' 
PRINT'ANCE FOR SOIL'* 

INPUT N1*S2*V2*S1 
S=S2*3. 14159/180 
V=V2*3. 14159/190 
MAT G=1DN(3*3) 

K3(0)=0 
FOR N=1 TO N1 

PRINT’ENTER NUMBER OF COMPONENTS IN* AND DEPTH OF LAYER* *N* 


1350 INPUT N2(N)fDl(N) 

1360 FOR M=1 TO N2(N) 

1370 PRINT'ENTER AREA-HORIZONTAL* NUMBER-HORIZONTAL* AREA-VERTICAL* * 

1380 PRINT'NUMBER-VERTICAL* REFLECTANCE* ANB TRANSMITTANCE FOR COM-* 

1390 PRINT'PONENT NUMBER' *M 

1400 INPUT M5(N*M) *H6(N'*M) *V5(N*M) *V6(N*M) *R9tN*M> *T<N*M) 

1410 J=2/3. 14159 

1420 A<N*M)=H5<N*M)YH6(N*M)>K(l-T(NtM) ) 

1430 A(NfM}=A(N*M>+V5(N*M)*V6CNfH)!t;(l-(R9(N*M)-H(N*M> >/2) 

1440 B(N>H) = H5(N*M)H!H6(NfM>*R9 (N*M ) FV5( N * M ) *VA < N*M ^ * ( <R9(N*M)-FT(N>M) )/2) 
1450 C ( N * M ) ^ H5 ( N * H ) # H 6 ( N ? M ) * R9 ( N * M ) 

1460 C ( N * M ) =C • N * M ) + J* VS ( ffl * M ) *V<b (N*M)*((R9(N*H>+T(N*M))/2)!«TAN(S) 

14/0 C 9 < N . M ) - J*V5 < N * M ) S- V6 ( N * M J I* ( ( R9 ( N * M ; +T ( N * M ) ) /2 ) *T AN ( S ) 

1480 C9 ( N * M ) =C9 < N * M > f H5 < N * M ) )KH6 ( N * M ) «T < N * M ) 

1490 K(N*M)-H5<N*M)i*H6(N*M>+JYVS(N*h!*V6(N*M)*TAN{S) 

1500 K2(N*M>=H5<NfM)4:H6<N*M>-FJ#V5(N*M'!t(V6(N*H)*TAN(V> 

1510 A1<N)-AKM»+A(N*M) 

1520 B1 (N)-B I (N)-t-B(N*M> 


1530 CKN.'^C I J-) FC(N*M) 
C2<N)-'"? ;• • ir?(N*M) 
1550 K3 (N; M(E • !K(N*M) 
1563 K4(N> !i'.2<N*M) 


OEIGINAL PAGE IS 
0P POOB QUALmr 
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:1570 NEXT H 

1580 K3(N)=K3<N-l)+K4(N)!HDltN) 

1590 K5(N)=EXP(K3CN-1>> 

1600 K6(N)=EXP(K3(N) > , , 

1610 G1<N)=SQR<A1(N)**2-B1(N)!S:*2) 

1620 Z1CN)=1/(A1CN>+G1(N)) 

IdSO Z2(N)=1/<A1<N>-G1(N) > 

1 640 Z7 < N ) =A 1 ( N > **2-K 1< N ) ##2-B 1 ( N ) **2 


1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1740 

1750 


IF Z7(N)=0 THEN 2500' 

Z 3 <N) = ((AKN)-K 1 (N) >ii:Cl<N)PBl(N)!J!C 2 (N>)/<Al<N)** 2 -Kl<N))t:* 2 -Bl(N):)c :»:25 
ZsInUKaiCnI+KICN) )*C 2 <N>+BKN)!«:C 1 (N) )/(AltN)** 2 -Kl(N)*)>! 2 -Bl<N)*# 2 > 

U1(N)=EXP<G1(N)*D1<N>) 

U2(N)=EXP(-lit:Gl<N)*Dl(N)) 

W3(N>=EXP(K1(N)*D1<N> > 

MAT E=ZER(3r3> 

E<lrl)=ZlCN) 

E<lr2>=Z2(N> 

E(1»3)=Z3(N) 


1760 E<2»1)=Z4(N) 


1770 E<2»2>=Z4(N> 
1780 E(2»3)=Z5<N> 
1790 E(3rl>=0 


1800 E( 3 » 2)=0 


1810 E(3r35=l 
1820 HAT D=ZER(3r3> 

1830 D(1f1)=UKN) 

1840 D<2»2)=Ur'(N) 

1850 D(3r3)=U3<N) 

1860 MAT P1=INV(E) 

1870 HAT P2=E*D 
1880 MAT P=P2>t:Pl 
1890 MAT H=G 
1900 MAT G=P!KH 
1910 NEXT N 

1920 A3=-1*(G(1 F3)-l*Sl)|t(G(2f 3)+G(3 f 3) ) ) 
1930 A4=G(lfl)-l*Sl^CG<2»l)+G<3»l) > 

1940 IF A4=0 THEN 2480 

1950 A0=A3/A4 

1960 R5(l)=0 

1970 R6(0)=0 

1980 I(lrl)=AO 

1990 K2t1)=0 

2000 I(3»l)=l 

2010 FOR N=1 TO N1 

2020 E(1t1)=Z1(N) 

2030 E(1f2)=Z2(N) 


2040 E(1»3)=Z3(N> 
2050 E(2f1)=Z4(N> 


2060 E< 2 f 2 )=Z 4 (N) 

2070 E< 2 f 3 )=Z 5 (N) 

2080 C< 3 f 1)=0 
2090 E( 3 f 2)=0 
2100 E( 3 f 3)=1 
2110 D( 1 f 1 >=UKN) 

2120 D( 2 f 2 )=U 2 (N> 

2130 n( 3 F 3 )=W 3 (N) 

2140 MAT P1=INV(E1 
2150 MAT P3=Pi*I 
2160 FOR M=1 TO N2<N) 

2170 Fl(NrH>=VS(NFM)H(06(NFM)!t:nT(NFH'+R9CNFMl )/2)>KJ^TAN<V) 
2180 Fl<NFH)=Fl<NFMmi 5 (NFH)l!H 6 <N 7 ii)*T' N fM) 

2190 F2<NfM J =M5 1 Nf M '< Si'76 f N f M ) < n i N fM J +R9 il N f M ) >/2 ) *TflN< 0 ) *J 
2200 F2(NfM) I'VNFHMHbiNFM'itiH/jlNFhl^RPCNFH; 

2210 F3(NFM>-<.<5<NFHAWf'(i‘FM:'!M I ! ( N f H ? f KV‘. M f M ) ' /2) 

2220 F 3 <NfM)=F I' !'AN ( 0 'iCTAW^ A) 
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2230 

2240 

2250 

2260 

2270 

2230 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 


F3 ( N » M ) =F 3 ( N r H > +HS <NfM>*H6(N»M)!t:R9<NTH) 

TUNfM)=Zl(N)*P3(lil >*(!/( K4(N)+Gl<N>))!*c<l-EXPaK4<N)+Gl<N) )*01<N>)> 

T2<NfM>=Z2<N)*P3(2f 1 )!):( 1/<K4<N)-G1(N) ) )*( 1-EXP << K4 < N >-Gl (N )) *D1 <N) ) > 
T3<NTM)=Z3<N)*P3(3f !>*<!/ (K1(N)+K4CN) ) )*< 1-EXP < <K1 <N>+K4 ( N) >#D1(N))) 


Rl<NTH>=Fl(NrM>*(Tl<N»M>+T2{NrM)+T3<NrH)) ^ . 

T4(N»M>=Z4(N))kP3(l r 1/!(;(1/<K4<N)+G1<N>) )#< 1-EXP( <K4<N)+G1 (N) )*D1(N> ) ) 
T5<N»M)=Z4(N)!KP3C2» 1 1/ (K4(N >-GKN) ) )sf:( 1-EXP< <K4<N)-Q1(N) >*DKN) ) ) 

T6<NrM)=Z5<N)iitp3<3T 1 )». ' 1/ (K1 <N)+K4<N ) ) )*< 1-EXP< (K1 (N>+K4<N) >*01 (N) > ) 

R2(N»li)=F2<NjM>*(T4<NTM>-i-T5(NTM)+T6<NrH) ) 
T7(N,M>=(1/(K1<N>+K4(N))>*(1-EXP<<K1(N>+K4(N))*D1(N)>> 


R3(N»H)=F3<NrM)!t;P3t3.1)*T7(N»M> 
R4 < N » M ) =R1 < N j M ) +R2 ( N » M > +R3 ( N r M ) 
R5(N)=R5<N>+R4<N»M> 

NEXT M 


R6(N>=R6(N-1)+K5(N)*R5(N) 


2380 HAT P4=E*D 
2390 MAT I=P4*P3 
2400 NEXT N 

2410 Q1=Z1(N1>*P3(1.1)*EXP(G1(N1)*D1(N1>) 

2420 Q2=Z2(N1>*P3(2j1>*EXP(-1*G1 (N1)*01(N1> ) 

2430 Q3=Z3(Nl>*P3<3rl>*EXP<Kl<Nl>*nl(Ni>) 

2440 R7=K6(N1 >*(ni+02+03) , 

2450 R=R7+R6<N1> 

2460 PRINT' THE RAIiIANCE IS’ ;R/3 . 14159? ’ 

2470 60 TO 2520 

2480 print’this f roblem has no solution because f rho(cll >-0. ’ 


wU i U •i-ta/aiv 

2500 PRINT'THIS IS THE REPEATED ROOTS CASE WHERE K=G , THIS SPECIAL* 
2510 PRINT*CASE WAS NOT CONSIDERED.’ 

2520 END 


READY 


OF SOOR 
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INF2 10J27 27-MAY-7S 


1000 PRINT "DO YOU WISH TO READ INSTRUCTIONS (YES OR NO)’J 
1010 INPUT A« 

1020 IF A$0“YES' THEN 1?30 ' 

1030 PRINT’ THIS PROGRAM DETERMINES THE RADIANCEf AS PROPOSED BY’ 

1040 PRINT ’GWYNN SUITS » FROM A SINGLE-LAYERED PLANT CANOPY WHOSE DEPTH’ 
1050 PRINT’ IS CONSIDERED AS INFINITE. THE PROGRAM WAS WRITTEN IN PARTIAL’ 
1060 PRINT’FULFILLMENT OF THE THESIS REOUIREMENTS OF JUAN H. CANTUj WHOSE’ 
1070 PRINT ’GRADUATE ADviSoR WAS DR, JOE CHANCE.’ 

1080 PRINT’ THROUGHOUT THE FOLLOWING CALCULATIONSr THE SAME UNITS OF’ 
1090 PRINT ’MEASUREMENTS SHOULD BE USED. THAT IS» IF THE DEPTH IS IN CM.f’ 
1100 PRINT’THEN THE AREA IS IN CM. **2. AND THE UOLUME IS IN CM.!t:*3.* 

1110 PRINT’ WITHIN A SINGLE LAYER. THERE CAN EXIST SEVERAL COMPONENTS’ 
1120 PRINT’ (LEAVES. STALKS. FRUITS. FLOWERS) EACH WITH ITS OWN SET OF* 

1130 PRINT’PHYSICAL AND OPTICAL PROPERTIES, THE AREA OF A COMPONENT IS’ 
1140 PRINT’DIVIDED INTO THE AREA PROJECTED ON A HORIZONTAL PLANE AND THE’ 
1150 PRINT ’AREA PROJECTED ON TWO ORTHOGONAL VERTICAL PLANES. IF THE AREA’ 
1160 PRINT’ IS MEASURED IN CM.*^2. THEN THE NUMBER OF HORIZONTAL AND VER-* 
1170 PRINT’TICAL COMPONENTS PER UNIT VOLUME SHOULD BE IN l/CH.*i(t3. THE’ 
1180 PRINT’TRANSMITTANCE (T) AND REFLECTANCE (R> OF A COMPONENT ARE UNIT-’ 
1190 PRINT’LESS AND 0<T<1» 0<R<1. 0<T+R<1. THE SUN ANGLE AND VIEW ANGLE’ 
1200 PRINT’ARE MEASURED IN DEGREES. BECAUSE OF THE CHOICE OF COORDINATE’ 
1210 PRINT’SYSTEM. THE DEPTH OF A LAYER SHOULD BE ENTERED AS" A NEGATIVE’ 
1220 PRINT ’NUMBER. • 

1230 PRINT ’ENTER THE WAVELENGTH IN NANOMETERS’! 

1240 INPUT J7 

1250 PRINT’ENTER SUN ANGLE. VIEW ANGLE. AND NUMBER OF COMPONENTS*. 

1260 INPUT S2.V2.N2 
1270 3=52*3,14159/180 
1280 V=V2*3. 14159/180 
1290 FOR M=1 TO N2 

1300 PRINT ’ENTER AREA-HORIZONTAL. NUMBER-HORIZONTAL, AREA-VERTICAL.’ 

1310 PRINT’NUMBER-VERTICAL. REFLECTANCE. AND TRANSMITTANCE FOR COM-’ 

1320 PRINT’PONENT NUMBER’ fM 

1330 INPUT H5(M) ,H6(H) .V5(M) .V6(M) .R9(M) ,T(M) 

1340 J=2/3, 14159 

1350 A(M)=H5(M)*H6(M)*(1-T(M) >+V5(M)*V6(M)*(i-(R9(M)-{-T<M) >/2) 

1360 B(M>=HS(M)*H6(M)*R9(M>+V5(M)*V6(H)*( (R9(M)+T(M) )/2) 

1370 C(M)=H5(M)*H6(M)*R9(M)+J*V5(M>*V6<M)*( (R9< M >+T (M ) > /2 ) *TAN( S ) 

1380 C9(M)=H5(M>*H6(M)*T(M)+J*V5(M))|:V6(M>*( (R9 ( M )+T ( M ) ) /2 ) *TAN( S) 

1390 K(M)=H5(M>*;!6(M)+J*VS(M)*V6(M)!fcTAN(S) 

1400 K2(M)=H5(M)*H6,lM>+J*V5(M)it!U6(M)*TAN(V) 

1410 A1(M>=A1(M-1)+A(M) 

1420 B1(M)=B1(M-1)+B(M) 

1430 C1(M)=C1(H-1>+C(M> 

1440 C2(M)=t.2(H-l)+C9(M) 

1450 K1(M)=K1(M-1)+K(M> 

1460 K3(M)=K3(M-1 )+K.2(M) 

1470 NEXT M 

1480 Zl=( (A1(N2)+K1:N2) >*B1 ( N2 ) *C2(N2 ) FBI ( N2 >**2*C 1 (N2) ) 

1490 Z2=B1(N2)**2+K:.(N2>**2-A1(N2)**2 
1500 IF Z2=0 THEN l680 
1510 Z=Z1/Z2 

1520 Gl=SQR(Al(N2)**2-Bl(N2>*)it2> 

1530 X1=1/(A1(N2)+G1) 

1540 X2=( (A1(N2)-K1(N2) ) *C1 ( N2 ) +C2(N2) *81 •' N2 ) > /( -1*72 ) 

15SC XJ-l/Bl(N2) 


■pA.G®_5 





1560 X4=< <AKN2)+K1 (N2) >*C2<N2)+B1<N2))CC1<N2) )/(-l!*:Z2) 

1570 FOR H=1 TO N2 

1580 Fl(H)=H5(H>!f;H6<H)i(tT<H)+V5(M)>l:V6(H)*(<T(M)+R9{H> >/2)W*TAN<V> 
1590 F2(H)=H5<H>*H6(M)*R9(M)+V5(H>!KV6(H)*< <T(H>+R9(H> >/2>*J*TAN<V> 
1600 F4<M>=V5(H)*06<H)*{(T(M)+R9(M))/2)!C<J»!»:2)»TAN(S>*TAN<U) 

1610 F3<M>=F4(M)+H5<H)*H6(H)NcR9(H) 

1620 Rl(M)=FKH>!t:(Xl*Z*(l/(61+K3(N2)))-!-X2*<i/<Kl(N2)+K3<N2))) ) 

1630 R2(H)=F2(M)*(X3*Z*<1/<G1+K3(N2) ) >+X4:it( 1/CKl <N2)+K3(N2> ) > ) 

1640 R3(H)=F3(M>»<1/<K1{N2)+K3(N2> ) ) 

1650 R4(M>=R4(H-1)+R1<H)+R2<M>+R3(M> 

1660 NEXT M 
1670 GO TO 1710 

1680 PRINT 'THIS IS THE REPEATED ROOTS CASE WHERE K=G. THIS CASE WAS 
1690 PRINT’NDT CONSIDERED. ’< 

1700 GO TO 1720 

1710 PRINT'THE RADIANCE IS" fR4<N2)/3,14159J ’ . ' 

1720 END 

READY 



;^OQ'R 
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27-HAY-75 


PRINT’DO YOU WISH TO READ INSTRUCTIONS (YES DR NO)‘r 
INPUT A$ 

IF A$'s.'^*Y£S* THEN 1240 

PRINT’ THIS PROGRAM DETERMINES THE RADIANCE FROM A PLANT CANOPY’ 
PRINT’CDNSISTING OF TWO LAYERS. THE DEPTH OF THE FIRST LAYER IS’ 
PRINT’FINITE WHILE THAT OF THE SECOND IS CONSIDERED AS INFINITE.’ 
PRINT’THE PROGRAM WAS WRITTEN IN PARTIAL FULFULLMENT OF THE THESIS’ 
PRINT’REQUIREMENTS OF JUAN M. CANTU » WHOSE GRADUATE ADVISOR WAS DR.’ 
PRINT ’JOE CHANCE.’ 

PRINT’ THROUGHOUT THE FOLLOWING CALCULATIONS* THE SAME UNITS OF 
PRINT’MEASUREMENTS SHOULD BE USCD. THAT IS* IF THE DEPTH IS IN CM.*’ 
PRINT’THEN THE AREA IS IN CM. **2* AND THE VOLUME IS IN CM.*Si3.’ 

PRINT’ WITHIN A SINGLE LAYER* THERE CAN EXIST SEVERAL COMPONENTS’ 
PRINT’ (LEAVES* STALKS. FRUITS* FLOWERS) EACH WITH ITS OWN SET OF’ 
PRINT’PHYSICAL AND OP'dCAL PROPERTIES. THE AREA OF A COMPONENT IS ’ 
PRINT’DIVIDED INTO THE AREA PROJECTED ON A HORIZONTAL PLANE AND THE 
PRINT’AREA PROJECTED ON TWO ORTHOGONAL VERTICAL PLANES. IF THE AREA 
PRINT’IS MEASURED IN CM.!K*2* THEN THE NUMBER OF HORIZONTAL AND VER-’ 
PRINT’TICAL COMPONENTS PER UNIT VOLUME SHOULD BE IN 1/CM.#*3. THE 
PRINT* TRANSMITTANCE (T) AND REFLECTANCE (R) OF A COMPONENT ARE UNIT-* 
PRINT ’LESS AND 0<T<1 * 0<R<1* 0<T+R<1. THE SUN ANGLE AND VIEW ANGLE 
PRINT'ARE MEASURED IN DEGREES. BECAUSE OF THE CHOICE OF COORDINATE’ 
PRINT'SYSTEM* THE DEPTH OF A LAYER SHOULD BE ENTERED AS A NEGATIVE’ 
PRINT ’NUMBER. ’ 

PRINT’ENTER WAVELENGTH IN NANOMETERS’} 

INPUT J7 

PRINT’ENTER SUN ANGLE AND VIEW ANGLE’* 

INPUT S2.V2 
S=S2)C3. 14159/180 
V=V2!i:3. 14159/180 
FOR N=1 TO 2 

PRINT ’FOR LAYER’ *N? ’ENTER NUMBER OF COMPONENTS’* 

INPUT NKN) 

IF N=1 THEN 1350 
GO TO 1370 

PRINT’ENTER DEPTH OF FIRST LAYER* * 

INPUT D1 

FDR H=1 TO NKN) . 

' PRINT’ENTER AREA-HORIZONTAL i NUMBER-HORIZONTAL* AREA-VERTICAL 
I PRINT’NUMBER-VERTICaL* REFLECTANCE* AND TRANSMITTANCE FOR COMPONENT* 

I PRINT’NUMBER’ f M 

I INPUT H5(N* H) f H6(N*M) * V5(N*H ) . VA(N*M > *R9(N*M ■ *T(N.M> 

» J=2/3. 14159 

• A ( N * M ) =V5 (N*M)>t:V6(N>M)*(l-(R9(NjM)+T(N*M))/2) 

> A(N»M)=A<N*M)+H5<N*M)it<H'i>(N*M)!C(l-T(N*H) ) 

I B(N*M)=V5(N*H)*V6(N*M)!lt( (R9(N*M)+T(N*M) )/2) 

> B(N>M)=B<NfM)+HS(NTM)*H6(N.M) *R9 ( N * M ) 

I C(NrM)=J*V5<N*M)*V6(N.H))(!C ( ( N * M ) +T ( N * M ) ) /2 > *TAN ( S ) 

) C(N.M)=C(N*M)+H5(N*M)*H6(N*H)#R9<N*H) 

) C9CN*M)=JTV5(N*M)^!V6(N.M))|t( < R9 ( N * M ) +T ( N . M > ) /2 > T T AN < S) 

) C9(N*M) = C9CN*M)+H5<N*M)!CHA(N.M)*1 (N*H> 

) K(NrM)=H5<N*n>*H6(N*M)+J!ityj(N.M>*VA(N*M)*TAN(B) 

) K2 ( N f M ) =H5 ( N * M ) *H6 < N * M > + J# V’J ( N . M ) VV6 C N * M ) *TAN < V ) 

) A1(N)=AKN)+A(N.M) 

) B1(N)=BKN)+B{N.M) 

) C1(N)=CKN)+C<N*M> 

) C ’CN> --E2<N) {-CV- N.M) 


OF POOR QUALm 
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1570 Kl<N)=Kl<N)fK<NiM> 
1530 K4<N)=K4<N)+K2«NrM) 


1590 NEXT M 

1600 K3<N)=EXF(K4<N-l)*ra ) 

1610 GKN)=SGIR<A1 

1620 Zl<N) = l/(Al<N>fGKN)) 

1630 Z2<N)=l/<AKN>-(jl(M> ) 

1640 Z7(N)=A.t<N>*>t'2-Kl'.N)!i:*2'Bi(N>**2 
1650 IF Z7(N)=0 THEN 2360 

1660 ^3<N)=<<A1<N)-K1(N))^C1<N)+B1(N>*C2(N))/(AKN>*#2 

1670 Z4(N)=1/BKN) 

1680 Z5(N) = < <AKN)+K1 (N) )>KC2(N)+Bi(N)!(<Cl <N) )/(Al (N)#<f2 
1690 HAT E=ZER(3,3) 

1700 E(1»1)=Z1(N) 


-Kl <N)**2-B1 (N>*#2) 
-Kl(Ni**2-Bl(N>*!»;2) 


1710 E<1»2)=Z2(N) 
1720 E(1»3)=Z3(N) 
1730 E<2>1)=Z4(N) 


1740 E(2»2)=24<N> 

1750 E(2t3)=Z5<N) 

1760 E(3fl)=0 

1770 E(3r2)=Q ' 

1780 E(3»3)=l 

1790 ON N GO TO 1800»1860 

1800 MAT D=ZER<3»3) 

1810 B(lrl)=EXP<Gl(l)#DU 
1820 D<2r2)=EXP(-l*Gl<l)i)!Dl) 

1830 D<3r3>=EXP<Kl(l)!KDl) 

1840 MAT Pl=E!t;D 

1850 GO TO 1870 i 

1860 MAT P2=IN0<E) 

1870 NEXT N 
1880 MAT P=P2Ji(Pl 

1390 IF P<2Fl)-P(2r2)=0 THEN 2350 

1900 XI (1)=(Z5< 1 >CB1 U)*P<2f 2)-P(2»3) )/<P<2» 1)-P<2»2) ) 
1910 X2<l>=~l^(Xl<l)+25a)*Bl<l)) 


19'^0 X3 < 1 ) =1 

1930 L1=X1 < 1 ))XEXP(G1 ( 1 ) *D1 )* ( 1/Bl ( 1 ) ) +X2< 1 )*EXP(-l!)!Gl ( 1 )!t:Bl )■*( 1/Bl Cl ) > 
1940 L2=EXP ( Kl < 1 ) #B1 > * < Z5 C 1 ) -Z5 (2 ) ) 

1950 X1C2)=B1C2)*(L1+L2) 

1960 X2C2)=0 

1970 X3(2>=EXP(Kl(l)!}cni) 

19BC FOR N=1 TO 2 
1990 FOR M=1 TO NICN) 

2000 FlCN»M)=05(NFM)W6(Nf H)*( (TCNrM)+R9(Nf M> )/27*J*TANCV) 

2010 Fl(N»M)-FlCNFM)+H5<NFM)*H6<N»M>)KT(NrM) 

2020 F2(NrH)=VS<NFM)^!06<N»M)!!!( ( T ( N t M >+R 9 ( N r M) >/2)*J!)!TAN(M) 

2030 F2CN»M)=F2(NrM)+H5<N.M)!CH6(NrM)*R9(N»M) 

2040 Tl<NrM)=05CNfH)*06(N»M)*( (TCN»M)+R9(N»M) >/2) 

2050 Tl(NFM)=Tl(N.M)!*t( J*:i:2)*TAN(0)5!:TAN(S) 

2060 F3<NrH>=H5(N»M)*H6(N»M)i*:R9(N»M)+Tl (NjM) 

2070 IF N=1 THEN 2090 
2080 GO TO 2210 

2090 01(N)=ZlCN)!t:XlCN)*(l/iK4<N)+Gl<N) ) )»(1-EXP< (K4(N)+G1(N) )S:D1> ) 

2100 Q2CN)=Z2(N)*X2(N)*( 1/CK4(N)-G1 (N) ) ))Jc( 1-EXP ( <K4<Ni-Gl <N) )*D1 > ) 

2110 CJ3(N)=Z3<N)*X3(N)*< 1/<K1<N)+K4(N> ) )!f:( i-EXPC (Kl CN)+K4CN> )*D1 ) ) 

2120 Rl<NrH)=Fl(NfM)*.<Rl(N)+Q2(N)+Q3(N) ) ^ 

2130 Q4<N)=Z4CN):CtXl<N)!(tC 1/(K4(N)+G1 (N) ) )*(1-EXP< <K4(N>+G1 <N) >sCBl ) > 

2140 Q5(N)=Z4CN)*X2(N)ii<<l/CK4(N)-Gl(N) ) )*( 1-EXP ( (K4(N)-G1(N) ) > 

2150 Q6CN)=Z5(N)*X3(N)^V l/C Kl ( N) +K4 ( N ) ) )C( 1-EXP ( (K1(N)+K4CN) >*D1 ) ) 

2160 R2<NFM)=F2(NrM)*(Q4(N)+0S(N)+a6(N)) ....... ^ ^ 

2170 R3(NrM)-F3(N»M>!CX3CN>.|t(l/(Kl(N)+K4(N) >)*<■ l-EXPC (l\l C N ) +K4 ( N > )*D1 ) ) 
2180 R4CNfM)=Rl(NFM)+R2(NFrt)+R3(NtM> 

2190 R5<N>-=R5(N)+R4(N»M) 

2200 GO TO 2300 

2210 01 (N)=Z1 <N '*X1 (N)>C(l/<K4fN)+Gl (N> ) ) 

2220 02(N)=Z3(N)*X3<N)!ir( l/(Kt <N)+K4CN> ) ) 
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2230 Rl(NrM)=Fl<N»H>iK(Gll<N>+D2(N>> 

2240 03<N)=Z4<N>!HXl<N)!i:a/(K4<N)+Gl<N>)) ’ 

2250 Q4<N)=ZSCN)i»:X3(N)*a/(Kl<N>+K4(N) ) ) 

22(40 R2<NtH)=F2<N»M)*<a3<N>+Q4(N) > 

2270 R3<N»M)=F3<N»H>*X3(N)!X(1/<K1<N)+K4(N))> 

2280 R4(N»H)=R1<N»H)+R2<N»M)+R3<N»H) 

2290 RS<N>=R5<N)+R4(NrH) 

2300 NEXT H 

2310 R6CN)=R6<N-1>+K3(N\)KR5(N) 

2320 NEXT N 

2330 PRINT'THE RADIANCE IS* »RA<2)/3. 14159f . 

2340 GO TO 2380 u_- ■ 

2350 print’this problem has no solution because h-x. 
2360 PRINT ‘THIS IS THE (REPEATED ROOTS CASE WHERE K=G. 
2370 PRINT'CASE WAS NOT CONSIDERED.' 

2380 END 


THIS SPECIAL* 


READY 
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